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ABSTRACT 


Voaveememtemor Une rav10 of the neutron-proton 


eect ric form factors, Ge /Gy, muere made from elastic 
ni 
Pee pIon-CecuLleron scattering to a precision of approx- 


imately 1 to 2% for the range of momentum transfers, a“, 


OryO.10 < qé < Og: me and for electron scattering angles 


muses) GG «120°, It was found that within experimental 
errors the RCAC eden ecnined from the ratio G, /Gp : 
agrees with the thermal neutron-electron Rete ar eae 
when relativistic corrections and proper deuteron wave 
functions are applied to the electron-deuteron scattering 


Results. 
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ie RODUCTION 


Ae EXPERIMENTAL OBJECTIVES 


iMicmemarce form factor Of the neutron, G as deter- 


ER? 
n 


mined by electron scattering is found to be essentially 


zero in the range of momentum transfers, ace OO res oS q* 


< 100 F-©. This fact is in apparent conflict with the 
results obtained by the scattering of thermal neutrons 

from atomic electrons. By methods to be described later, 
measurements of the slope of the neutron charge form factor 


at very low momentum transfers of q® ~ 1077 nee show that 





dG, | 

n = +0.0193 + 0.0004 F¢ Giza 

a(q®) | 2 
q 70 
It is very doubtful that a non-zero slope for Gr at q°+0 
Pe 
ean be reconciled with G, = 0 at q® > 0.3 F7*. 
nN 
Drickey and Hand [1] in 1962 measured Gr /Gy for 


i 


7 to an accuracy of about 2%. 


values of q*, Ones q° ege He 


Their results yielded Gp = 0, but when reanalyzed by Casper 
n 
and Gross [2] (1966) using recently developed deuteron 


LoraIbivEsuoLe C@rrecbLions and new deuteron wave functions, 


the value of G was brought into fair agreement with the 


i 
n 


value of (1-1). Thus the entire problem is further compli- 
cated by the extreme sensitivity of Gp to the adifferent 

n 
deuteron models assumed. Hence, if very accurate measure- 


Men tse OL Ge. Or q° small, can be made by electron-scattering 
n 


el 





merema1ques, a comparison with the value (1-1) could be made, 
which might yield an answer to the "correct" deuteron model 
emoe chosen. 

It is the purpose of this experiment to (a) demonstrate 
the basic agreement between the thermal neutron-electron 
slope value of Gp and the electron scattering Gr 


n n 
(b) try to determine which deuteron wave function (model), 


pect KG 


best describes the deuteron. 


B. NUCLEAR STRUCTURE FROM ELECTRON SCATTERING 

Pence une Carly 1950's production of beams of high 
energy electrons by linear accelerators (LINACS) has made 
Eemeolole the study of the electromagnetic structure of the 
mueweon. It is expected that the dimensions of any particle 
Oeeiiass m is of the order of n/mc, which for a nucleon is 
about 2 x 10724 em. Hence, to investigate the structure in 
this size range it is necessary to use a probe having about 
these wave lengths or shorter. lLinacs produce electrons 
with just such wave lengths. The 100 MeV NPS Linac thus 
"sees" only the outer structure or surface of the nucleus. 

Plleectron scattering Nas several definite advantages 
over other methods of nuclear structure research. Briefly, 
the important reasons that electron scattering is such a 
useful tool are the following: 

(a) The basic interaction between the electron and 


target nucleus is known; the electron interacts with the 


charge-current density and magnetization density of the 


2 





Meteloe oince the electromagnetic interaction is rela- 
tively weak, of the order of a = 1/137, the fine structure 
constant, measurements on the target nucleus can be made 
taenouy greatly disturbing its structure. 

ites in COnlrast CO the situation with protons, 
Meweronis and alpha particles, since they are strongly inter- 
acting particles, and the scattering mechanism cannot be 
Miweye separaved from structure effects in the target. 

(bd) Photons interact electromagnetically with the 
Miers, but Since a real photon has rest mass zero, it can 
Meaiorer Only a detinite amount of momentum to the struck 
Mmoereus.. This 1S in contrast to Secor scattering where 
the momentum transferred to the nucleus is limited only by 
Mremciicrey lost to the nucleus. 

itempeenaoistLy that an incident electron of total 
energy E is scattered elastically at an angle @ (i.e., o, 
the cross section for elastic scattering at angle 6) can 
be calculated precisely on the assumption that the target 
nucleus consists of point nucleons. That means the charge 
and magnetization densities are delta functions. The 
actual cross sections differ from the predicted cross 
section if the charge and magnetic moment are Spread Oulice 
amempie ciiterence is usually in the direction of smaller 
values of the cross section since the smeared out structure _- 
fremiceauehe field against which the particles scatter. The 


deviation becomes greater as the parameter, q, the momentum 


es 





: 


transferred to the struck nucleon increases. The (three) 


mementum transfer is defined by 


+ + aS 
eee = Da ae. 


where De and Pp aromuncmrnec| Gene, and final values of the 

electron momentum. The effect of structure of the nucleons 
memiormally expressed by the use of form factors, G, which 
meemoeliecvyedq to be functions of the structure and momentum 


meanster only; 


aH ey = SQoscmsecelOuelOrsseavtering trom physical particies 
a cross section for scattering from point particles 


(1-3) 
Note that similar ideas occur in all scattering situations 
Mamviiech the Suructure of the scattering center has an 
milmence, €.f., X-ray scattering form factors reflecting 


the electronic structure of atoms. 


Peete seTRON SCATTERING FROM DEUTERIUM 

The deuteron, as the only bound two-nucleon state, has 
Peonmpne source Of much information for nuclear physics, 
due initially to the great interest in the two-nucleon 
problem. Originally the two-body problem was thought to 
be of great importance because of its relative mathematical 
simplicity, and because of the hope that the inter-nucleon 
force was additive, so that knowledge of this system's 
solution would allow extension to more complex nuclei. It 
is now generally assumed that the nuclear force is not 


B@@ative (esee,ethe nuclear force, unlike the electromagnetic 
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force, "saturates") and hence the two-nucleon interaction is 
Merenecessarily the best source through which understanding 
Sec me nature of the nuclear force can be gained. 

Since free neutrons are unstable and decay with a half- 
ieee Ot |/7 minutes to a proton and electron, deuterium is 
Pmeolvenlent source for neutrons as targets. HO~. cAas 
experiment deuterium in the form of solid deuterized poly- 
ethylene (CD,) 5 was used as a target. Assuming that 
deuteron wave functions are known from theory, neutron 
structure is inferred through a comparison of electron- 
meswemeand €leceron-deuteron scattering results. Conversely, 
knowing the nucleon form factors from theory makes compari- 
sons of different deuteron wave functions possible. 

There are several possible outcomes to an electron- 
deuteron scattering experiment: 

(1) e+ Dw. et + D! (elastic scattering) 

(2) e+ D-7+e!t +P+N (inelastic scattering, break-up) 

(3) Radiative scattering 

(4) Reactions with pion production. 

Now (2) occurs only above 2.2 MeV excitation energy and is 
not of interest here; (3) can be calculated; (4) is not a 
problem for incident energies smaller than 150 MeV. The 
theory for process (1) was developed initially by Jankus 
fi (i'>5o)°in the first Born approximation, which should 
be adequate because of the small deuteron charge. The 
deuteron is described by nonrelativistic wave functions 
while the electron is treated relativistically. Jankus's 


result will be presented in Section II, part B3. 


US 





PELempta ate a relativistic deuteron theory were initi- 
ated by Jones and carried forward by several authors, prin- 
Smeclly Gourdin, Adler and Gross [4, 5, 6, 7, 8, 9, 10, 11). 
Meouwrhus OF their attempts will be presented in Section II, 


peamee B55. 


D. THE NEUTRON-ELECTRON INTERACTION 

The neutron-electron (n-e) interaction is the term used 
to describe part of the electromagnetic interaction between 
Mgemiacurron and Che electron. This interaction may not only 
Hema specific interaction between neutron and electron, but 
Mem~cen any Charged particle and the neutron. Furthermore, 
there exists the possibility of some n-e interaction of a 
nonelectromagnetic (weak) nature. For the present case 
however, results that stem from the interaction of an elec- 
tron with the internal electromagnetic structure of the == 
neutron will be termed "the neutron-electron interaction." Cee 

' Since Miemmeluron has Zero total charge there is no 

-Coulomb interaction between a neutron and an electron. 
Both particles have magnetic moments and hence there exists 
a spin- dependent magnetic dipole - dipole interaction 
between them and also a velocity-dependent interaction 
between the magnetic moment of the neutron and tie Magnevire 
iiemidmeasseeciated with the convection current of an electron 
imenovulonwe ihnese interactions have been well documented 
and are not of interest here [12]. If there are regions of 
non-zero charge-density in the neutron, (i.e., a charge 


Soperatrvomeexlata)) then any charged particle probing the 
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neutron will experience spin and velocity- independent 
PeeecuroOstatic forces, These are the forces of interest in 
the n-e interaction. 

At least some charge separation in the neutron is 
expecved on the basis that the free neutron possesses an 
mmemalous magnetic moment. This contribution to the inter- 
eermon 1S called the Sere Cte or Foldy term [13]. The 
remainder of the interaction (due to charge separation) is 
Wgemresult Of the fact that a neutron can be part of the 


enema lssociated into a negative pion and a proton, 
n<«> p’ + 1 


If the neutron spends 80% of the time as a neutron and 20% 
of the time as a proton + negative pion, then an electron 
penetrating the 'neutron' would see an electric field 
sieemetvn equal to that produced by a charge of 0.2e, e the 
Droton charge. The force resulting from the field would 
have a very short range since the tm would screen the pos- 
itive charge at distances greater than the 7 p* separation. 
(The separation can be on the order of the pion's Compton 
wave length h/m_C x 1071 em.) The force should be 
attractive and can be represented as an extremely narrow 
potential well. The purpose of all the experiments 
described below was to determine the well depth Vo? escen— 
tially by measuring the scattering amplitude. 

Wrewfirst, of three techniques to measure the n-e inter- 


action was introduced by Havens, Rabi, and Rainwater [14] 
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momero4ty. It consisted of scattering low energy neutrons 
(~O.1 eV) from liquid lead and bismuth, noting that the 
meee SCavvering cross section of an atom includes an 
observable term arising from interference between the 
Mileweus and the electrons scattering coherently, while at 
higher neutron energies (~10 eV) this interference term 


is almost absent. Results on liquid bismuth yield 


dG. /a(q*) = 0.0225 + 0.0070 F*. 


n aa) 


The work of Hughes, Harvey, Genders. and Stafne [15] 
depends on the fact that the difference in the refractive 
indices for neutrons at an interface between bismuth and 
liquid oxygen (liquid mirror) comes largely from the neutron- 

electron scattering amplitude. The measurement of the angle 
mee eel exvtcernal reflection gives the difference between 

the two scattering amplitudes and hence the n-e interaction. 
Here the results are 


aG. /a(q°) = 0.0200 + 0.0019 F°. 
E 2 
mA eg. 0 


The most recent and most accurate measurements have 
Soewmobwerohn anc Ringo [16] (1966) utilizing a technique 
Oterermi and Marshall [17] (1947). It depends on the fact 
MiotiaescaLveering Neutrons, the interference between the 
scattering from the nucleus and from the electrons in an 
atom leads to an asymmetric angular distribution of 


scattered neutrons in the center-of-mass system, This 
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mecurs because the atomic form factor depends on the scat- 


mering angle. 


A summary of the Krohn and Ringo results are 


smown in Table I. 


Gas 


Argon 
hevoton 
Xenon 


Weighted 
mean 


Vi fev) 


-~3780 + 250 
-~3800 + 130 
~3670 + 90 


-3720 + 90 


Adio: a 


dG, /a(a*) Fe 
nN 


eo = 000 12 
CeO Or 0007 
Oe O19 Ort 20 0005 


0.0193+ 0.0004 


otatistical 
weight 


Os 
0250 


0.57 


Momiceimoorvant to note that all measurements agree that 


dG. /a(a°) is not zero or negative. See Appendix A for 


Nn 


dG 


m@iemeerart's of the calculation of Bn from ee 


a(q*) 


1, 





ieee Lu bORY 


A. THE NONRELATIVISTIC DEUTERON 
moe Ground State of the Deuteron 
There are three possible states of the two-nucleon 
system, the di-proton (pp), the di-neutron (nn), and the 
deuteron (np), of which only the deuteron is known to be 
Stable. Experimentally measured static properties are: 
(1) Binding Energy W = 2.226 + 0.003 MeV 


(2) Magnetic Moment Hq = 07057356280 .00002 nm 
where lnm = 1 nuclear magneton = ee and 


m,, is the proton mass 
(3) Electric Quadrupole Moment, Q 4 


= 28.2 ee where 1F = 107? em. 


DEC oe ore 


(4) Spin Angular Momentum, S = 1h 


(5) Radius gee = 4.3 F (see page 23 for definition of Pq) 
(6) Mass Mz = 2.01410 umu (unified mass units) 
= 9,50 Fol 


Because the neutron has zero charge, the nuclear force 
binding the deuteron cannot be electrical and neither can 
it be gravitational, for this force is far too weak to 
pmevide a 2.2 MeV binding energy. Tighe. aunueled Iz thoree 
Piece ts Serong and attractive is postulated. By ignoring 
some experimentally determined parameters a simple quantum 
mechanical treatment of the deuteron is possible. 

Assume that the nuclear force is central, i.e. the 
interaction potential between the neutron and proton is 


some function V(r), where r is the separation of the two 
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nucleons. This assumption is in conflict with (3) above 

because a non-zero quadrupole moment implies a nonsymmetric 
Siete Gisuribuvaon. Nevertheless, assume the central | 
potential, and take the ground state of the deuteron to be 
ae: Hence the Schroedinger equation in the center-of-mass 


system is, 


_ne 2 => > > 
con w(r) + Vir) plr) = Ev(r) (2-1) 


where the reduced mass is, 


m= —P D2 2 C2) 


and Mm,» m, are the masses of the proton and neutron. The 
average of these masses is M, and r is the relative vector 
coordinate between nucleons. E is the energy of the rela- 
tive motion of the two particles. Since the ground state 


is spherically symmetric (S state), w(r) depends only on 


=> 
mer). so that 


eon (223) 


ie 


w(r) 








The qguantity u(r) is called the reduced radial wave function. 


Now using (2-3) and (2-2), (2-1) can be rewritten as 


2 
a Meee ue) = Eucr) (2-1) 
d“u(r) 
where u(r) = <= and & = -wW = -2.226 MeV. For the 
ar 


Square—well potential 


Za 





(2-4) becomes 


Case L r< Py 


Mier) focutr) = 0 # £(2-6a) 


where 


M M 
a= {—(V_+E)= | 3(V_-W) 
he 0 fn? 0 (2-6b) 


The general solution of (2-6a) 
is 


u(r)=A sin ar+B cos ar 
(2-7) 
Requiring that wW(r) be finite 


eber>O0 and wW(r) = 0 at r>»” 
fesecs B = 0 in (2-7), hence, 


aS) 


moe) = wesin ar 


Now at r=r 
O 
Pous (2-0) and (2-11) state 


Ce e 


A sin ar 
O 


and 


-Yr 


-yCe : 


Agcos or 
O 


both w and w' (the 


(2-5) 
Case. li r > Le 
u"(r) - y“u(r) = 0 (2-9a) 
where 
Y = a (2-9b) 
he 


The general solution of 
(2-9a) is 


Lt) Con ope. 2a 0) 


the boundary condition at ~ 
requires D = 0 so that 


u(r) = Ce Yt (224) 


slope) must be continuous, 


(2 ade.) 


(2-13) 


Mme racio of (2-13) for (2-12) yields 


MW 
z 


a cot aor_=Y7ye= - 
: h 


ee 


(2-14) 





fie is a transcendental equation relating the range of the 
potential r, to its depth ioe: Typical solutions of (2-14) 


Or r, vs depth are given in Table II. 


WABI 








Depth of 


Range r. (F) Potential (MeV) 


1.0 120 
iS 5 
2.0 36 
2a 25 
eS Ses 


The results of this rough calculation are quite good. 
Independent measurements of ry yield numbers on the order 
met Fermi, If ear e.0 F is accepted, the size of the 
deuteron can be specified by 1/y = 4.3F, which is about 
twice that of the range a On thes poucnbial. This explains 
Why the deuteron is a loosely bound system (its binding 
energy is 1.113 MeV/nucleon, vs. 8.0 MeV/nucleon in the 
average nucleus). 

There is no experimental eviderice for any excited state 
of the deuteron. It can be shown by simple theoretical 
mesumenvs that in order to produce a bound excited state 
the potential well depth must be much greater than in the 
Pound stare. It must be concluded that no bound state 


exists for 2 # 0 in the deuteron. 


eo 





e. Review of Low Energy Neutron-Proton Scattering 


Certain parameters pertaining to the deuteron have 
Sect established by investigation of the unbound n-p state, 
Pm@eo 1S through neutron-proton scattering experiments. It 
memclear that these parameters should be predicted by any 
Beecrcneory Of the deuteron. With this in mind this brief 
review is given. 

The n-p scattering cross section has been examined 
closely, both at high and low energies. The cross section 
Pivem@sesourongily on the energy of the incident neutrons, 
ev eetiereies, below 10 MeV, the scattering is essentially 
due to neutrons with angular momentum 2 = 0 (S-wave scatter- 
ing). Thus the angular distribution of the scattered 
Mamerons 15 150Cropic in the center-of-mass system. Since 
been provon has molecular binding energy of about 0.1 eV, 
only neutrons with incident energy greater than 1.0 eV are 
eemsrderca.,. Henee the neutrons strike essentially free 
Preevons. 


The Schroedinger equation for the scattering problem is 


Vew(e) + s LE -v(F)| H(P) = 0. 
fi 


feta central povential and for neutrons incident along the 


Z-axis the above becomes 


aceon oo [E- v(r)| Ceca eee (2-15) 
ae 


waich has the asymptotic solution 


ay 





ikr 


lea, e 
€ 





V(r ,&) A + £(@) 





(2-16) 


r+ © 


The f(@) is called the scattering amplitude representing 
the scattered neutrons (waves); ng represents the incom- 


one plane wave. It is not difficult to show that 


oT lece) |, (2-17) 


Using separation of variables, w(r,@) = R(r) Y(@), (2-15) 
G@ivides into a radial and angular equation. The angular 


Dart is 
atk & (sine ax | Puce) YY =O C28) 


with solutions Y 9 (@) = P, (cos @), the Legendre polynomials 


momomder &, since the potential is central these solutions 


are independent of the potential V(r). The radial equation 
is 
wl (r) + pe - u(r) - EY J uicr) = 0 aie 
r 
where 
Ke = a and U(r) = emu te) 
fh n 


Pek r= uiry/r. Let the potential V(r) be "turned off." 


fimen the solutions of (2-19) are again the spherical Bessel 


functions eget) Tt can be shown that as r > » 
. Simi kr tae) r 
j, (kr) er ae os Voted [ 


a) 





and hence 





Uy (r) sSin(kr - &1/2). (2-20) 


ra @ 


fiemey is “turned on" but r is still greater than rr, the 


solution of (2-19) must still correspond to that of a free 
particle, and hence be the same as (2-20) except for 


possibly a phase factor. Thus 


‘ QT " " 
Ol) saa Sin kr mi + by) Weron . 2p) 


ties eeneral solution to (2-19) for "V on" at r > ie is 


sin(kr - £1/2 + ©2) 


v(r,0) = > Cy - 


C 


[M1 8 


P, (cos®) 
g 


0 
(2-22) 


» B,P,fcos&). 
g=0 Ce 


But this must be the same as (2-16). Now using a Rayleigh 
expansion on ge in (2-16), and writing all quantities as 


complex exponentials, equating (2-16) and (2-22) yields 


16 
Cc. = qe eine A 
eee eS 
: 
anda 
i 16 > 
r(e) = L) (2n+1)e * sin&,P,(cosé). (2-23) 


R 


The square of the scattering amplitude is the differential 


scattering cross section (2-17), hence 





_ f ao(e) .. . ay inc peo) 
Om -f qa sin? dé 2 (2£+1)sin Oo. ( ) 
Q 
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Note that S-wave scattering (2 = 0) is predominant for 
em@ersies less than 10 MeV. By taking the range of nucleon 
force to be 2F and assuming a 0.5 MeV incident neutron 


(in center-of-mass frame or 1 MeV in the lab) it can be 


seen by forming the ratio of the coefficients B, in (2-22) 
ena t 
Demic 
i uD °C 
a Wisse) 
O 








This states that only about 9% of the scattering is due to 
feterons with £ = 1. 
In (2-24) let 2 = 0, then 


e } Z 
fee 2 _ eo 
Om = me) Sol on = i(— ° (2-25) 


Perk = O, sin 6G approaches Oo: imomder thas Om remain 
finite 6 /k must approach some constant, called -a by 


convention. Thus 


On = bra C2=26)) 


foe rowa 1S called the scattering length. Consideration of 


eiae Solution 


6 
: ae ee: meee - 

u(r) ~ sin(kr ~ 6) = Seana tee) =a area!) 
demonstrates the fact that the asymptotic solution becomes 
vl omnoelimesand that the scattering length, a, is the 
micercept On the r—axis and is obtained by extrapolating 


the radial wave function inside the potential, well beyond 


ay, 





the range of force ae (oUuCewmmaliustraves the sienifi- 


mimee Of the sign of the scattering length. The SCauver ime 


Pength a is positive if the state is a bound state. 


(a) Figure 1 (b) 


Positive scattering length 
indicates a bound state 


Negative scattering length 
indicates an unbound state 


The energy dependence of the scattering cross section 
ean be expressed in terms of the scattering length a, and 
another parameter which has the dimensions of length and 

is known as the effective range. Denote this parameter 
by p, where 0 is of the same magnitude as the range of 

nuclear force oe Errective range theory predicts the 
phase shifts 55 as a function of energy. Consider S-wave 
scattering only. The wave equations for two energies 


a Hoke /2m erate) O are 
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Comin 0 Case II V = 0 
ut (k® - SY yy = 0 (2-27a) vit 4+ Ky = 0 (2-294) 
N : 

u" - eae 0 (2-27b) v" = 0 (2-29b ) 

O H eB O 
mietiply (27a) by ug, (27b) by Meets eee ay sb Vice 2.91D)) 
Meand subtract, getting by Veand subtracy, setring 
a ees eo ac ' eee 
AD (uus - uu )= k uu, (2-28) Bee Nee Ve) = k aw (2-30) 


Subtract and integrate the 


Gur Susp’ - vv' t+ vv) 
O O O O 


Because u = v, u. = v_ for 


Oo O 
jeft hand side contributes 


Furthermore u(0) = u, (0) = 


! ? 
(vv. SR ) 


r=0 





In analogy with (2-20) the 


u(r) v(r) =C 





Oe 


where C is a normalization 


1 - 
meas cot O. = 


above equations. Then 





r=o D = 
=k (uu_-vv_)dr. (2-31) 
O O 
r=0 
0 
eae ro? Lhe expression on the 


HOMimt at ume Upper Jimmie. 


0, hence (2-31) is 


2 | Cillie—) 7 acre . 
O O 
0 


solutions u and v can be written 


ae) 


sin(krté.) 
Seeal@aai On) = Sa 
O 
constant. Now (2-32) becomes 


oo 


2 
k / (uuo-vv, )dr 
@ 


and in the limit of zero energy neutrons u = u 


>? a O 
that 
2 
_-l k 
k cot 6 = = EO a 


a 





where 


oo 


_ 2 2 
ae = 2 | (vo = u,)ar 
0 


ms the effective range. 


Hence the cross section is given by 


Mor 


On = rere (2-33) 
k Fa Gao) 
Besides being a function of Ke = omE/* , the cross section 


is expressed in terms of the two parameters, effective range 
p and scattering length a. The cross section is independent 
of the.form and shape of the potential. Both po and a can 

be determined by measuring the cross section at different 
low energies of the incident neutron. 

The theoretical estimate of the value of the scattering 
eross section for &£ = 0 is far below the experimentally 
Observed value. In fact for a square well of reasonable 
width and depth 


On * hiS x 10724 em* 


while 


O . = PA oe 1072" om*. 
experimental 


This conflict was removed when Wigner realized that the 


theory for o, is based on a potential whose "strength" is 


- 
calculated from the deuteron, where the spin is 1, i.e. 
neutron and proton are parallel, while scattering experi- 


ments are done using randomly oriented nuclei. The 
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ee 





scattering experiments then represent an average over 
Perallel] (triplet) and anti-parallel (singlet) states. The 
Weight factor for the singlet and triplet cases are 1/4 and 
3/4 respectively, hence the total cross section for scatter- 


ing is 


ae 1 
2" 7% + F %s 
where O, and Oo, are the triplet and singlet cross sections. 


Using (2-33) and working backwards assuming o, = 20.4 barns, 


C 


ra, = 4.3 F, and ee el eereeon be determined that 


ya = 0.066 MeV with Ae = 25 F, where the scattering length 
a. is negative. Hence the singlet state is not a bound 


deuteron state. 


poeeemeonomenological Two-Nucleon Potentials 


On the evidence presented so far it can be concluded 


5 
S4 state 


perce the solutions to the Schroedinger equation for a 


that: (1) The deuteron is predominantly in a 


central potential are quite good; (2) since Hs # He + Hs 
and Qy x 0, there is reason to guess that noncentral forces 
are acting (the inequality of magnetic moments in a complex 
problem and could well have other interpretations); (3) the 
forces between neutron and proton are different in che 

mere let and triplet states, implying that nuclear forces 
me Spin—dependent. There is a point of confusion on (3). 
Spin-dependent does not imply that this contribution is not 


@®@envural in nature. ittemeans thaw the usual central scalar 


merential is different according to the relative orientation 


Se) 





mops, 2 Situation that has no direct analogy either in 
meecuromagnetism or atomic physics. 

A realistic deuteron potential must incorporate these 
ideas and more. The potential must be symmetric in the 
coordinates of the two nucleons and invariant under rota- 
means and reflection of their space eee ainates. This 
means the potential is a scalar. 

Some of the more detailed features of the nuclear 
potential can be directly related to the spin and isospin 
dependence of the nucleon interaction. An important 
Semrrtoutrion for the nuclear symmetry potential comes from 
mae exchange character of the nuclear forces. Thus a 


serber exchange potential, Bae fOor=thnew central interaction 


would be 
¥, = = (1 +P") V(r) or 
yo = + 1 - + (1+ 7,-4,)(1+6.°6,)1V(r) (ney) 
c 7 3 i 1°72 1% 


ae -> -> > : a 2 - 
where Tis To and O45 % are the isospin and Pauli spin 
matrices respectively. 
The occurrence of a rather strong spin-orbit force in 


mee nuclear interaction gives rise to a spin-orbit coupling 


mothe two body force; 


ats 


jl > > = > (2 7 
Vig = = Vile) (ry - 25) % (Py - Po) (Sy +55). (2-35) 


fis Vis contribution is based on Wigner's idea that 
velocity-dependent forces are acceptable when they depend 


only on the relative momentum p of the two nucleons. 


Se 





All these facts taken together yield a general potential 


VE VoG) a Vi) 6-6, + Vi(r)S,,+V,g(r)l:8 (2-36) 


Paere Bae mepresenus the —Spin-spin interaction and Sy5 is 
mae vensor operator 


— 


ee ce, ©) = (c,-c5). (2-37) 


The form of the tensor operator can be inferred by analogy 
with the classical problem of the potential energy of an 


electric dipole d. in the field of another dipole ds at a 


1 
distance r.* Here 


Si -3 = a ae < a ae 
Vag 2 304, r)(d,+P) - a, d. |. (2-38) 


The r7> comes from the fact that the electric charges 
Mizeract with a potential proportional to ie But the 
r-dependence of "mesic charges" which are thought to be the 
source of nuclear fields is unknown. Hence the r3 depen- 
dence must be replaced by an arbitrary function of r, Vin(r). 
inemstmueture of the nuclear force is closely related 
to the properties and interactions of the whole family of 
aerouely inveracting particles. It has not been possible 
Seetar to derive the forces between these particles, or 
their masses, from simple assumptions regarding the basic 


structure of the strong interaction. Nevertheless, certain 


relations between the interactions and masses of the 


*See PopehiGabassei Or properties Of the Lensor operator 


S10: 
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strongly interacting particles can be established. Of 
particular significance is the relationship between the 
nuclear force at large distances and the pion-nucleon 
moeeraction. 

The force at large distances can be described in terms 
of an exchange of mesons, in a similar manner as the elec- 
tromagnetic interaction can be analyzed in terms of photon 
exchange. The special role of pion exchange is due to the 
m@eriness Of the pion mass as compared with that of other 
strongly interacting particles. The interaction associated 
with the exchange of a particle of mass m_ is limited to 
a range of the order of the compton wave length A, =h/m_c. 
In fact, the intermediate states involved in such an 
exchange have an energy of at least m_ of and so are limited 
Gena duration of order h/m, 0°. During this time, the 
emitted particle cannot travel farther than h/m,c and thus. 
the interaction is expected to decrease strongly for 
distances greater than dae If n particles are exchanged 
samultaneously, the corresponding intermediate energies 
free y) X m, and the range is dfn. (ie winveraceson abe) tne 
Pmeesoucdistances is thus determined by the exchange of 
Single mesons. 

The asymptotic form of the one-pion exchange potential 
(OPEP) is uniquely given by the mass and symmetry properties 


of the pion (J’ = 0°, T = 1) 
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=A 


eo (4. t (5, +3 Seon pee 
Vopep = 3 ne Mao §7y' To) [65 FD) 2 : (2-39) 
where x = m_ c/h = 0.70 poe r</he Orvis sO OO 2 : 


TT 
E@ewstrength of the potential is determined by the coupling 
menstanct f for the process N+ N +m. The meson field is 
pseudoscalar since g" = (ee aidmisevecror since Lf = I. 
Equation (2-39) reduces to (2-36) for the deuteron. 

Early definitive work on the phenomenological nucleon- 
nucleon potential in the OPEP approximation was carried 
out by Glendenning and Kramer [18]. They attempted to 
construct triplet (spin) - even (parity) potentials that 
Poems movOcle to the OPEP and are modified in the inner 
regions so that the deuteron properties are obtained. The 
results of an 8-parameter fit to an equation similar to 
(2-39) are found in Table III, entries 1-9. Coefficients 
mee cet OOLvenyial are varied such that the empirical values 
of Qy = 28.2 F© and ae = 5.44 + 0.02 F are reproduced. A 
hard core radius ry was chosen so that the potential has a 
bound state at W = 2.226 MeV. (The hard core is discussed on 
page 38). 

Many other authors have contributed to the data [19, 20, 
21]. One of the most useful has been the Hamada-Johnston (HJ) 


[22] potential given by 


> > 
V= Vor) of Vin(r)S,5 a Vigh’s + Vip fbi (2-40) 


5D 





TAS IG ES iL 


Moti fication r,(f) a, (F™) p(F) Q4(F*) Py th) ; 
GK1 Poo fomebearky es 26.79 7.415 0.02714 
2 0.2466 5.368 1.709 28.21 6.281 0.02693 
3 ie oebomminole 828.17 5.970 0.02595 
4 ie peso oten (O03 5920.22° 7.103 0.02677 
5 inccocmoccomeiy705 26.02 6,710 0.02663 
6 0.3924 5.384 1.726 28.12 5.957 0.02676 
7 iP oscomeomigfel. colo. 2e.14 6.029 0.02564 
8# ie O0o Senin yo0 20.18 5.622 0.02654 
g* iie2 oes sein i> 26.29) 7.425 ° 0.02670 
HJ 0.485 - ieee oso 6. oY 6 0.02656 
BL 35 - - - 4.60 - 
Experiment : isn | = - 


mers the asymptotic D state/S state ratio 


* GK8 and 9 include the spin-orbit term in (2-34) while GK1-7 
mee without it. 


+ N. XK, Glendenning and G. Kramer, Phys. Rev. 126, 2159, 1962. 


Tt G. P. Auffray, Phys. Rev. Letters 6, 120, 1961. 
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where 


tH 
| 

Q 

at 


2 


> > > 
(Coeac ec) ain) 


and contains central or tensor Vin» spin-orbit V and 


line 7 


second-order spin-orbit V COUPONSMGS ett 1s from this 


LL 
mevential that the later used Partovi wavefunctions have 
memecerived, The radial functions are restricted by the 
condition that at large distance the central and tensor 
meecntials should be described by the OPEP. The radial 


forms are: 


Ve = vo(%, °F) (85°55) YC) [1+ a ¥(x) +b Y2(x)] 


Vin ss ea OAC, ) Z(x) [1 +ap¥(x) + dp¥*(x)] 


. O 
Vig o Ere Voz (x) [1 +b). ax) | 
V = 9 V WX) 44 Y(x) 4 65 ¥“(x)] 
LL ic O xo LL tole 
pen oo 
ve are es = 3,65 MeV 
mc -x 


Xx =— re eo ne. = 
h 1.43F ° x 


In addition the potential has been assumed to have a com- 


ponent giving rise to an infinite repulsion at radius 


fee = 0.343 co ee Eee eo 


Si) 





The values of tne quantities ae Vin» etc. have been 
memectllacved for moeur different types of potentials. The 


results are shown in Table IV. 


TABLE IV 





V (MeV) Vin (MeV ) Vy g (MeV) Vip (Mev) 


Singlet even -~1460 _ - aih\2 
Triplet even —- 207 -642 34 668 
Singlet odd 2a - - -6688 
Triplet odd - 23 173 -1570 -1087 


Values of HJ potential at r, = 0.485 F. 


The "hard core" mentioned above is a phenomena that 
Cannot be investigated through low energy nucleon-nucleon 
scattering. In the neighborhood near and above 200 MeV 
repulsive components in the nucleon-nucleon interaction are 
seen. One way of introducing a repulsive interaction (which 
might account for nucleon saturation) is in terms of a very 
strong repulsive potential of very short range - "the hard 
core." The infinite short ranged repulsion is only one, 
rather extreme, way of accounting for the observed effects. 
Pmive repulsive potentials, "soft cores," as well as inter- 
moions depending explicitly on relative vélocity can account 


for some experimental results. 
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A fairly recent idea by Feshbach and Lomon (FL) [23] in 
1967, is to replace the interior region by appropriate 
boundary conditions at some finite radius. As long as the 
interactions responsible for these boundary conditions 
involve energies much greater than the bombarding energy, 
it may be expected that the boundary conditions will be 
approximately energy independent. The FL model also employs 
Phe (TPEP) two pion exchange potential. 

4, Solutions of the Non-Central Force Deuteron Problem 

solutions to the non-central force deuteron problem 


are obtained by substituting the potential V ince the 


OPEP 
schroedinger equation and solving for the wave functions. 
Angular momentum and parity arguments can determine the 
deuteron ground state angular wave function. 

The total angular momentum of the deuteron is J = 1 
made it 1S in a definite parity state. The value J = 1 can 


be obtained from different combinations of the orbital 


Seerar momentum &, and the spin s such as: 











et ash 
s = 0 aa 
ol , 
s = 1———"P, 
R= 2 s = ] on 
gmqgd no others. If the deuteron is a mixture of these it 
can only be s, + om Or oer + 2B since the combination 


muse Nave a definite parity. Since an almost spherically 
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Seammerric ~round state exists, the obvious choice is the 
Meetbatre combination. It is not difficult to show that the 
odd parity states Gey a oDeo) yield the wrong values for 
Q4 and Hg. 

The angular momentum wave function Y eae for a two 
particle system in the orbital angular momentum state 2 and 
mmm state ©, forming the total angular momentum J, with 


Z-component M is given by 


ae <2smoM_JIM> vy,” "s (2-41) 
4 WSJ > TR's g Xs - 
MM 


M M 
where the Yo are the spherical harmonics, ce eet deus Dall 


wave functions and <& Sm, M | JM> are the Clebsch-Gordon 
coefficients. Two spin 4 particles have the following spin 
ewaves: 


Triplet (s = 1) 


x7 elec Iya (2) na 1 ja ( 2) (2-422) 
Xq = <''s74%s|10>0(1)B(2) + <tts4s7¥s| 10>0(2)8 (1) 
= =[a(1)8(2) + a(2) 8(2)] (2-42b) 
v2 
xp ts <tas7eo|1-1>8(1)B(2) = B(1)B(2) (2-H2c) 


Singlet (S = 0) 


[ (18 (2) : B(1)a(2)] ei) 


M(H 
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where a = ie and g = i) are "spin up" (proton) and "spin 
down" (neutron) states of the nucleon. The arguments of 
the a's and 6's stand for particle number one or two. The 
angular momentum wave pees) ants for the two states 


Mader consideration are 


i 
a <0101]11> Yea(1)0(2) = Yo(1)a(2) (2-44) 
011 0 0 
De jek, Mes 19) So) ees 
1-—M_ OM 
Pa -) <21(1-M,)M_]11> ¥ : one 
S 
Boke Soe ok. 
= 2 v8 acayecay /B v3 Fala (18(2)+8(2)4062)| 
iL le 
+ fa Y5 a(lja(e2). (2-45) 


The Schroedinger equation in the center-of-mass system 








is 
=i a es 6 = Ed (2-46) 
om C ABSA ye 
where 
= eer) il w(r) 1 
eo Dp or You +r Y o43 (2547) 


is the ground state deuteron wave function: For ease the 
linear spin-orbit and quadratic spin-orbit parts of the 


Meeential (2-40) have been ignored. Rearranging (2-46) 
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2 
Cs v4 v, - EYP Y onl 011 OY | 
n Ve : 55; + el Yor : 24 13] OnE 2= 10) 


In writing (2-48) the following properties of S have been 








12 
used: 
$194 011 2 E43 (2-49a) 
1 
and S10 Yard = v8 as : Ba (2-49b ) 
Note that 
eee 2) Le (2250) 
oe ar’ ~ 722 
where 
2 9 
L* = h. sin ¢ a (sin ¢ Se) + 5; (2-51) 
sin'¢@ og 
and 
m m 
L°Y,” (6,9) = £(2+1)n° ¥,” (6,6). (2252) 
Hence 
vey d = bey 2 zee) 
LL r O11 
ar 
ana 
y2@y 1) 2 (1 aw _ 1 2(eer) , (2-54) 
r f211 a2 Yor: 2 


Ho 





+ 
Multiply (2-48) vy Y a integrate over the solid angle, 


M 


mmo make use of the orthonormality property of theY pss 


imo. £eCT 


z a 
p u"(r) + (V, - E)u(r) = -¥8 Vy w(r). (2-55) 


. 
Similarly, multiply (2-48) by Y ae and integrate to get 


a 
om 





|w" (xr) - 6 we) + LV Wee E | w(r) 
r 


= . /8 Vpu(r). (2-56) 


These coupled differential equations were first obtained by 
Rarita and Schwinger [24]. 

A method for solving these equations has been described 
by Hamada [21], who states that (2-55) and (2-56) have two 
sets of linearly independent solutions. Thus choose (Cu, Wy) 
and (uy ,W5) with the asymptotic requirements that the first 
Solution behave as an S wave function and the other as a D 


wave function: 


f eer 
ih (2-57a) 
Wy QO 
and Us 0 | , 
2-57b 
We 2 el 1 + 3/ar + 3/(ar)*| 


Starting with the behavior at large r, integrate (2-48) 
inward to get the two solutions for arbitrary r. Then u(r) 
and w(r) are obtained as a linear combination of these two 


solutions: 
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u(r) 1 u, (rv) uy (r) 


SI 5 +h : (2-58) 
w(r) Cho) = w, (r) w,(r) 
The requirement that u and w in (2-58) vanish at the core 
radius te gives an equation for Po which is solved by an 


interpolation procedure. The asymptotic D state/S state 


ratio A is then determined from 


u,(r,) w,(r,) 
Aro ae 


The solutions (2-58) are found in practice by numerical 
methods. Now having the wave functions in hand the magnetic 
Moment and quadrupole moment can be predicted by taking the 


expectation values 


<u> -[o un dr, (2-60) 

| 

ees gt Or ease (2-61) 
and 

<Q,> =] >" Q, oar, (2-62) 

Q4 = aa a ro (2-63) 


ieaeecal results are listed in Table III for various poten- 
tials. (See Appendix C for example calculations of re 


and <Q5>-) 


B. THE SCATTERING OF ELECTRONS FROM PROTONS AND DEUTERONS 
herwuuena li lyeecleectron Sscacvtvering cross sections have 


Beem handled through the techniques of first order 


yy 





time-dependent perturbation theory, the main result of which 
is the Fermi Golden Rule No. 2. This rule relates the 
mes ition probability per unit time, [, to the square of 


Mmeeeiiaurix Cléememt for the transition. Specifically, 


2 ea 
ea! | | 41> | 


aN 
dE (2-64) 
where GN/dE is the density of states like |f> around the 


energy E, and Hy PeeleminnMerachromeananl Itonian, The cross 


section is then Gefined as 


. 


~ UIneident particle flux’ (2-65) 


Quantum electrodynamics (QED), which is the interaction of 
MeemOirac spinor field with the electromagnetic field, pro- 
vides the mechanism for computing the cross section via 
S-matrix theory. In general terms the S-matrix is a unitary 
operator that transforms the initial state |i> into the 


fapel state |f>, i.e., S|i> = |f>, where 


oo 
yee Ae > = cA? falixy.. false Be) 28, 
ake (2-66) 
eet 1S the chrenological operator that orders the inter- 
action with respect to causality. The expression (2-66) is 
integrated over all four space and hence is not physically 
measurable. If Sey is normalized to unit space-time volume 
iem@eit can be written 
Sri] | 


VT 5° (Pp - Py) I Mey | 





(2267) 
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where 5" foeowleoumwaimeisteonal Gelta function expressing the 
conservation of four momentum in the transition from |i> to 
Mm, The cross section for a transition from a two-particle 
Mittal state |i> into some group of final states is in the 


center-of-mass 











2 
mea 27 ) il 1 S [ 4 2 > > 
Soh 7) 2,1 6 poy) Mpg" GP] + AP 
p ECS.) E(2.) final final 
Py Po? soins momenta (2-68) 


meee rration over suitable variables yields the desired dif- 
ferential cross section do/dan. 


At the heart of (2-68) is the matrix element |M He 


pil: 
Rueeme mine Characuer of the interaction is found while the 
other factors represent kinematical quantities. The form 

of M es can be surmised by examination of the Feynman diagrams 
mere che process under consideration. 

Much of QED is conveniently described by Feynman diagrams. 
Here the diagrams will be mainly used to supplement the 
Pooerstanding of the electromagnetic interaction, but their 
memes CorPeainiy not limited to this role. If the diagram 
memoc Craw, then there exist specific rules for construct- 
ing the S-matrix therefrom. 

i -omeOlwsceuverine 1s sraphically represented by 
Figure 2, Both graphs are equivalent, but (a) is most often 


feem since if is just a little more suggestive of the 


Scattering process in the center-of-mass frame. 
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(a) (b) 
Figure 2 


The time axis is directed upwards unless otherwise noted. 
tee NOorizontal axis can represent a space or momentum inter- 
meme GQistinction between the two being of no importance 
mee. lhe arrows on the lines point in the direction of 
ieereasing time for a particle. If the arrows were reversed 
the particle then represented would be the antiparticle of 
the original one. Figure 3 sometimes is said to represent 
particles in the "scattering channel" while the "crossed or 
annihilation channel" is represented by the same diagram 


rotated by 1/2. 


rotated ccw 


by 
T 
Zz 
scattering annihilation 
channel channel 
Figure 3 
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Lines beginning and ending within a graph represent "virtual 
particles." Table V below lists some of the more important 


Meaph elements seen. 


TABLE V 


Graph Mactor in Lhe 
Element S-Matrix 


Leo a eae) annihilation (creation) 
E PJs 1 Pp of an electron 
p Dp 
en los annihilation (creation) 
Bo Ce BD p of an antielectron 
(k) 
ra ai annihilation (creation) 
20, Ou pam loOron 


Interaction 





eu 
—_ Sein Culeer ron. Profaror 
p +m 
ae LY NS JL 
— PiOLon proltavor 
q 
1 
ae ae meson progator 
eet an) 


Y 6p !-k) photon-electron vertex 





phovon-nucleon vertex. 
hip Or meross slab ened 
k yFors (p—-p!—k) circle represents ignor- 
agee or the true anver— 
action; mathematically 
this ignorance is repre- 
SecmpComeyeLormnm Lactors 
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Paeeelascvlcer oe uron—-) rovon ocattering from Point Protons 

To calculate the elastic electron-proton (e-p) scat- 

tering cross section from point (Dirac) protons in the meet 
Bern approximation it is assumed that (1) the electron is a 
physical point, and (2) the one photon exchange process des- 


cribes elastic scattering. The graph of Figure 2 is then valid. 





Figure 4 


feeerice figure define the kinematical quantities P5oPp> 


PsoPp as the four-momentum of the incident and final 
reerronl and proton respectively; q is then the four- 
momentum transfer 


oe 


°C 
q = (ps a Pp) 


= (P, - Palee (2-69) 


Choose the convention that ae < 0 for electron scattering.* 


By looking at Figure 4 it can be seen that 


*™e metric used is a* = q* - 4°. In the lab frame the 


recoil energy given the struck particle is small_with respect 
to the total incident electron energy thus qe = maieue Oe aie 
four vector q¢ is then said to be spacelike. In the annih- 
ime. On channel qe iz 4Me_ > O and is timelike. Another 
metric commonly used is ee ote Tiecmcanhereny metrics 
yield cross sections and form factors with some sign dif- 
ferences. This possible ambiguity is being remedied by use 
of the variable t for +q¢ in many newer papers. 


Ho 





eae = ee (2570) 


Nolte 


which represents the interaction of the electron four- 
current Jy with the four-current of the proton g’ via 
maeecval photon progator ao. 

Equation (2-70) is the equivalent of taking the inter- 


action Hamiltonian to be 
Eeea = iene (x) (2-71) 
where a¥(x) is the electromagnetic field of the proton, 


Ween 1n turn depends on the proton current TeGe- Regard- 


less of the point of view the cross section (2-68) becomes 


e 


O 
y eS u 
o ~ 6 (pp- Pp, t Pe- Pid, aad | (2272) 


Smmoredicts the form of the current for the electron, 


J 





U(Pp)Yulp,) (2-73) 


and by analogy for the proton 


Sie 





Oey UCP, (2-74) 


where u,U and u,U are the incident and final electron and 
p@eren Spinors, respectively; and Y is a 4x4 matrix. Sub- 
aeeeition Of the currents into (2-72) and carrying out the 
Pal followed by a sum over final 


spins and an integration over final momenta will finally 


trace calculation of |[M 


preld the result 
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++ 
o\e ox -1 ne 
do re COpmae 2 ae ina 2 8 
(5) es ee (u a Sin : (- - tan _ 
wirac 5 


Sin MC om 
Erovon 
(2-75) 
where 
=) SOG = do 

=) ———- = Mott cross section = | — (2-76) 

(i sin! & (Sie 
= E (EIS 


-] 
and = sin® 2) = recoil factor 


are the results of the kinematical quantities in (2-68). 
The square of the invariant matrix element IMp.| yields doe 


eer ac proton 


Boe 
a 1 - 24> tan* 3 


a (225) 
om e 


Mes 
In obtaining (2-75) the approximation i 0 since m,< E 
has been used. Equation (2-75) is correct for scattering 
from a point proton with no anomalous magnetic moment. 
But the proton is known to have an anomalous moment so 


the structure of M must be changed to reflect this fact. 


ft 
imacdi tion, the possibility of nucleon structure should 
bewctaken into aceount in Mes (the electron remains a point). 
Figure 4 is redrawn to schematically illustrate the change 


a iat Mea 


Sel 





rEReure ss 


The bubble at the photon-nucleon vertex represents the non- 
locality of the photon-nucleon interaction and the anomalous 


moment of the proton. The matrix element is now 





; ost 
Dia Ot, Y 
q 
where 
U — u 
J U(P,.)0 UCP, ) (2-79) 
and 
yp _ ee 2 ON, 
Or = F,(q yy" + F(q NG) Gar (2-80) 
ine quantities Fo and B, are undetermined real functions of 


a’, while ob =e amie: een Foldy [25] has shown that while 
(2-80) is not unique it has the proper general form. 


With the above changes (2-75) becomes 


| Ce 
a z =) rex [Po- GS fa(ry +e P)* tan® Sere 
R Mott 4m,¢ Pp p 
(Beco) 
where the final bracket is the invariant matrix element 


nee, 1°. and where ne 1.79 is the protons anomalous (Pauli) 


Sy 





moment. Equation (2-81) is known as the Rosenbluth cross 
section. The factor F, (a) PmEouna ,O represent the 
mmemon's Charge Gistribution and its normal magnetic moment, 
while F.(q°) represents the anomalous part of the magnetic 
Moment. Fo and Fy are called the Dirac and Pauli form 


factors respectively. Yennie [26], and later Hand, Miller 








Paeewitson [27] pointed out that Fy and Fe, should not be 
femstaered to be electric or magnetic form factors. More 
meeperliy certain linear combinations of Fy and Py defined 
to be 
Ac 2 
= hoo 1 (2-82) 
E il Nm-e* Diee ae 
and 
Cee = Het Kal ap i eal (2-83) 
M a} De a q> + 0 p p 
are termed Sleeumcnana marnetice form factors. The major 


advantage of (2-82) and (2-83) is that in the Rosenbluth 


cross section (2-81) the cross term (Fy + KPa) is removed, 


9 Ae ep ty 
Cie ai a > 
io a (<2) ae. = a. 2 tan* = ae (2-84) 
Mott ie fa. Mme 
ume 


Pam Ouimerie —ecOsogesecl1on formulas the form of the 


é =. It can be 


quantity in braces is A(a®) + B(q°) tan 
shown that this is the most general form for the one photon 


B-wellamee process, regardless of the structure of the 


2 





Zmeaerirome and tie Pprovon. The interpretation of A and B in 
terms of Gp and Ga Oc Fo and Fy Gewoenas Ol COuUrSe on the 
eeeucture assumed, 

eee Nucleon Form Factors 

In the text above no real attempt has been made to 
Generently explain the form factor concept. So far they 
meres Deen characterized mainly as correction factors that 
emu itiplied’®y The scattering cross section for point 
(structureless) particles yield the cross section for 
meeeectineg tromaphysical particles. Thus, as introduced, 
these dimensionless quantities are meaningless unless some 
meee cLULONS canibe made from a basic theory of matter about 
Pfemseructure and behavior of the form factors. 

While the simple idea of the form factor representing 
some physical extension in space can be illustrated by non- 
aol vistic Cxamples, a more comprehensive discussion of 
their meaning needs the help of the so-called dispersion 
memaptons. Both viewpoints are investigated below. 

In Section II, part A2, the scattering amplitude f(@) 
fewer ouna 19 a pertial wave solution of the Schroedinger 
“Howlonm 1Or the scattering of neutrons from protons. Of 
interest here is the solution of the Schroedinger equation 
for the scattering of electrons from two potentials V(r), 
Shewwloh a point charge source, the other with a smeared 
Siitechargse adistrabution. 

Vi G(@eahimdsagumooLons as im section II part Ad, the 


Pcm@rPocodimeer eguetion for scattering is again 
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(Vo + ko) ¥(r) = Bv(r)v(r). (2-85) 
A 


mero differential Schroedinger equation is solved formally 
homey making (2-85) into an integral equation. The inverse 
eee ae linear differential operator ce F K~) roe aneintlegral, 


hence 


Vi(r) -[ c¢r,2") 3 Veco Vir’ jidr" (2-86) 
f 


where r' is a dummy variable of integration and G(r,r') is 
met i1om Characteristic of the form of the differential 
operator. G(r,r') is known as a Green Function. The 
Beneral solution to (2-85) is the sum of the solution to 


the homogenous problem (in this case it would be the 

> > 

unscattered plane wave ee since oP v(r) = 0) and some 
h 

Parercular solution to the inhomogenous problem which is 


given by (2-86). Hence 


V(r) = a + fore.) 2 Vir')V¥(r')ar'. (2-87) 


Hoes ts Called the integral equation for scattering, and 
meemmeOoluLion is found by iteration. The first iterated form 


momealled the Born or plane-wave approximation 


> > 
ai erent 2m 


» {G(r,r') elk Pr apr, (2-88) 
a 


Comparing (2-88) and (2-16) yields the scattering ampli- 


mide * 


*In order that the comparison be made the correct form 
of G(r,r') must be determined. 


Dp 





> > 
mee 25 fet-* V(r')dr' (BaRo 
A 


where q =k. - Kp. For a spherically symmetric potential 


Pee (2-89) is 


f(@) 32, | iu seenmete svicre') dirt, (2-90) 
h'q 0 


ferme Point charge of Ze, the potential V is zZe“/r, thus 
2 
f (6) = mo (2-91) 


while for a distributed charge of Ze of the form 


Zec -ar 
ie >; ‘l-e Dacia = Os 
Da 8 2 
ome a a 
f.(e) = — =. => =f (¢) —-— . (2-92) 
D noa? cana ite. q? ae 


Hence the structure of the distributed charge shows up as a 
mioworm Limes the point charge scattering amplitude. It can 
Pememown Chat the form factor here is the Fourier transform 
of the charge distribution. Figure 6 illustrates the 
behavior of (2-91) and (2-92) as a function of q°. The 
slope of f(@) at q* = 0 is proportional to the radius of 
the structure under consideration. Before continuing it 
must be stressed that all that has been said so far of form 
Boerors reters to the non-relativistic case. In a rela- 
tivistic example difficulties are immediately apparent. A 


maeere distribution, spherical in one reference frame is 


elongated in another frame so that the meaning of the radius 
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ee slope = O= structure 
; radius 10 


slope # O= finite 
structure 
f(¢) radius 


Figure 6 


Seeeeecharge distribution is not clear. No longer can a form 
maewor be Chought of properly as the Fourier transform of 
Mmeemetotributvion, One has to be concerned with the form 
fevers themselves and not their three-dimensional Fourier 
transforms in a particular reference frame. 

the Meson theory of nuclear forces mentioned briefly in 
Pewerens | part D and TT part A3 states that virtual pions 
Meemeonwocvantily being emitted and reabsorbed by a nucleon. 
fueemer, two nucleons are stuck together through an exchange 
eemerous. "Now, Since nucleons are surrounded by a cloud of 
penis, the virtual photon emitted by an electron scattered 
Peomeune nucleon might have in fact, interacted with the 
Diemeconm meson cloud. So the simple electromagnetic inter- 
action of a photon and a nucleon may be quite a bit more 
complex and depend on the dynamics of the strong interaction 
Seec@e plon and the nucleon. Dispersion relations then 


attempt to understand form factor behavior as a function of 


oa 





G imeverms Ol the Strong interaction dynamics. What 


follows is a simple minded attempt to make plausible the 
mer Of the equations put forth Boe wie “form Lacwor , 

The basic idea of dispersion relations (spectral repre- 
Sentations, or S-matrix theory of strong interactions) is 
tCiae a quantum mechanical amplitude for a physical process 
is the boundary value of an analytic function of a complex 
variable. For electron scattering F(q°) would be the 
boundary value along the real line of the analytic function 
F(z) where z > qs + i¢. A dispersion relation is an inte- 
gral equation that relates a dispersive process to an 
absorptive process. The best known example is that of the 
dispersion of light in a dielectric, where the complex index 


Of refraction n(w) = n_,(w) + i a(w) is expressed as an 


ae 
2.) 
integral over all frequencies involving the absorption 
coefficient a(w): 


CO CO 


ae w'Im n(w) _ a(w')dw' (2-93) 
Ce = 2 12 a) | te 2 
9 - Ww 9 = 0 


30 


(P represents the Cauchy principal value of the integral) 


A classical optical example of (2-93) is Sellmiers equation 


for the anomalous dispersion of light, 


: oe 
a = i= > ae (2-94) 
A” -dA 
i=1 i 
: ; : la 
where A is the incident light wave length, d. the i 
resonant wavelength, A. acensvanc., and N the number of 


natural frequencies of the medium. 


58 





in electron scattering 


2. al ImP(a*) 
F(q') = = Some Co (2-95) 
hme q 
TT 
is a direct analog of (2-93). Now three mesons, the 0, $¢, 


and the w-meson have been found experimentally To have the 
Correct properties so that they may couple to the nucleon 
aembO Che Virtual photon. They correspond to the N natural 
frequencies above. Hence it is plausible to write keeping 


(2-94) in mina 


3 
Som 
Fae) = 5 5 jeer (2-96) 
i=] Lo = 
me 
SL 


In terms of the physical form factors G, (2-96) might be 











written 
‘es Swp Sop Spp ‘ 
= + + t+1l-g -g, -g 
a. 2 2 2 wp “op “pp 
p coe ee ae 
m< m< m< 
u) d 0 
and 
2, ya a ee 
ccd) = aa ae 5 + 1- Uy Uy ca ® 
: =: oa 
MA m4 uP 


Memes tne @'s and W's are coupling constants to be fitted to 
experimental data. The twelve coupling constants in (2-97) 
are not all independent and their number can be significantly 
reduced. The best fit available for the proton was used 


here to predict G, for use in Gao meee Vries [26] adjusted 
p 


0 





the free parameters appearing in (2-97) so that a minimum 
momootained in a yo aemcOmcnmcriamenval data. The followtng 
values were determined: ep OIC Eg, SOS a | 
aap 1.191, where m,° =a 6 wes, my” = 26.6 F-* are 
experimentally confirmed resonances of fairly narrow width. 
Since the p resonance appears as a rather broad peak (750 +t 
100) MeV and since some observers [29] have shown that the 
effective position of the po mass can be at about 600 MeV, 
the mass of the p resonance has been taken as a free 
parameter, determined to be m,* oon eee toremcan= 
megane it should be noted that both (2-94) and (2-97) are 
approximations of (2-93) and (2-95) respectively in the 
sense that the resonances considered are taken to be delta 
Pumevions in the integrands of the more accurate form factor 
expression (2-93) and (2-95) 

It is interesting to note that the theory outlined above 
actually preceded the discovery of the mesons that it 
meauitres. Dispersion relations as applied to the strong 
interaction was developed in 1957 by Chew [30] et al on the 
Deasis Of a two-pion exchange approximation. This theory 
fPeomeu able Lo explain form factor behavior, and it later 
showed (1959) that a three pion resonance must exist in 
order that theory and experiment agree. 

iiimiicmoy Ort ermehactor notations exist today. 


Coupled with their often misleading names and differently 


Gefined static limits, they present a baffling array. 
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Pommonly, une Svarving place is with the two functions 
F,(q°) and F,(q°), called the Dirac and Pauli form factors 
(really there are four form factors, two for the proton, | 
two for the neutron). They describe the structure of the 
meme Dirac coupling and the anomalous Pauli coupling 
between the nucleon and the electromagnetic field. Bo and 
Fy Pieniiou directly related to the electric and magnetic 
Wemenmc distributions of the nucleon, but in older papers 
mieverare Often erroneously called electric and magnetic 
form factors. Still another general name for Pec Fy 
emeweroex form factors referring obviously to their use 


meee Current operator ge. The Dirac and Pauli form fac- 


meres nave the following static limits: 


ee) 
a A 


ll 
ke 
es 
o~ 
© 
Ww 

il 
ae 


F, (0) = 1.79 = « aco 
p n 


sae (2-98) 


Gaon MNibler and Wilson [2/7] popularized the use of the 

physical form factors G. They were introduced to simplify 
PaemnOsenbluth cross section by removing the cross product 
terms in F, and Fy and are in no way more fundamental than 


the vertex form factors. They are defined 


ae oe CH 
Gy = Fy + KF = (1tk)Fyac (2-99) 
Eee 2 2 

where n = -q /4m", m the neutron or proton mass, kK the 


eiomalous proton or neutron moment. Bon and nag were used 
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briefly before the G factors became popular. The form 
ma1ecvOrs Gy and G_, do merit the names, electric and magnetic 


form Tactor, and they have static limits 


HH 
© 


Gp (0) = 1 Gi (0) 
p | n 


I 
Hl 


Cue) = 2.792 Caen] 1.0] (2-100) 


m m 
p 1G) 


Peers Oven esound that still another type of form factor 
memeoonvenient to work with. In isotopic spin notation 
where eS i Cmetne proton and lege -~l1 for the neutron, the 


Peapee and magnetic form factors can be written succinctly 


morepoth the neutron and proton as 


_l us 
= iment) Gy, eee) Gy, 


p nN 
oe: i 
Gi Ms ail + eo Gy ae sy Gs 
p n 
er as, 
ae Ae 
Gp = 5(G, +G, ) + 57,(G, -G, ) 
p n p te 
1 oa 
G ==(G. +G )+=y7y /(G@ -G_ ) 
m 2 wi m De az: m,, Mm. 








1 
G" = =(G + G. ) = 
E 2 = ED q® , 2 
S 1 ih 
@> = Se = (Gee) =(y +p) Cz VOl) 
m EZ M, Mm. q@ maa 2B n 


and the isovector form factors as, 
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ae 1 
Ge = —(G. = G_ ) = 
EB 2°", E ie, ° 
ay elte -g ) os (2-102) 
m ee m 
p n Ca 0 


Mese (iso) form factors are convenient for use in theoreti- 


@meework Decause the electromagnetic current i may be split 


S V 


mmo tsoscalar and isovector parts, l.e. J. = J + ley 


| Hu Hu 
The intermediate states in the photon-nucleon vertex can 
be characterized by this isoscalar-isovector idea in that 
the two pion resonances with (T=0, J=1) w and 9 manifest 
Ss 


Peremoeclvyes in the isoscalar form factors Ca and Gs , 


whereas the third resonance (T=1, J=1)p contributes to the 


Tsovector form factors G. and oe ites Inveres tance To 
S 


note that only the isoscalar form factor eGh = (Gh + Gy ) 
appears in the deuteron cross section, which might a ae 
Premeiiace Lie ceuLeron current is an isoscalar. 

There is an important relationship between the physical 
form factors that moeeriscovered ime olewat Orsay by 
feamenn, Taylor, and Wilson [31]. In low ac measurements 


(actually the work was at q® = 2,98 poe) of the e-p cross 


section the following was noted: 


G m 
aD (2-103) 


Mater data verified this up to approximately 100 F~°, and 


(2-103) came to be known and called the "scaling rule." 


There is some suggestion of (2-103) when the vertex form 
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BeewOroware exXPressea in terms of the physical form factors, 


for example inverting (2-99) yields 


iL a a ]-7 


The above have poles at n=l unless (2-103) is true. Drickey 
and Hand have experimentally verified (2-103) to include 
the deuteron with 

G M 


i = 1.71 = a Ug: C2=1045 


G G 

m Me 

pe GA, (2-105) 
Vo Un 


Figure 7 shows that the scaling rules can be used as an 
established empirical fact (to at least 20 a) without 


Paerocucing appreciable error. 


eG. 
2.79 Gry 
i@| 
os ta i;y; 
0.5 
Oerao a1 0 20 30 5 
Bae . ORSAY 
foam ord 


a4 Harvard 


ede ile 
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Teelecblcen Lectaon Weutberon Scattering 


Following Enereane general procedure as Rosenbluth 
did with the proton, Jankus (1956) derived [3] the cross 
section for elastic electron-deuteron (e-d) scattering from 
memos. deuteron in the first Born approximation. His 
feomwlts ditfer qualitatively from the proton cross section 
MicmrinoomLans Way, 1.e€., there are three form factors 
Miverecuerizing the deuteron while two suffice for the proton. 
This is in agreement with the work of Glaser and Jaksic [32]. 
ier study showed that the cross section for scattering a 
Per vistic electron from a potential with spin J, contains 
wee torm factors. The three deuteron form factors are 
found to correspond to the charge, quadrupole, and magnetic 
ements in the static limit ac See ie seeneral forme! tne 


fees Section is as it should be for the one PMO rome <change 


2 6 
: 


peeccitacally Jankus found that 


process, Ata") + B(q*) tan 





an AQ 2 Q 


-l 
ae (22) 2 a “ pea? + = n“D<(q°) 
Mott We 


d 
+ én { 1+2(12n) tan* 4 od?) (2-106) 
where 
D,(q) -| [ue (r) + we(r)|5o(4)ar —1 (2-107) 


0 


fmemene charce SiLructure ractor, 


65 





2 2 
Da (4°) = “3 nf fuCr)w(r) - | j(t)dr —~ 50 (2-108) 
0 


fmomeune Quadrupole moment structure factor. 


co 


Dy (a) -{ [ue(r) - gw (r)]5j(rar + + | [u(r)w(r) 
0 0 


4 w°(r) 








: ]dp(t)ar—-1 - 3P, (2-109a) 
3 


mecounus for the contribution of the intrinsic magnetic 
moments Of the proton and the neutron to the scattering 


Brocess and 


oo 
eee, CS @ : a: 
Dy(q") = af Ww (r)[5,(1) feo t))| dr + Ph (2-109b ) 
0 
fmameme Maenetic contribution to the scattering process 
oeeene from the convection of charge in the deuteron. 


Together (2-109a,b) are 


one M E | 7 
Dy(a°) = >| + W)2Dy + Dy (Silo 
miso Biven are the static limits (a° +0). The ie and Jo are 
epnerical Bessel functions, T = lal r/2, Qo, tne, Quadrupole 


moment and Py the percentage D-state (see page 127). 

While the currents in e-p scattering were easily found, 
the deuteron current presents somewhat more of a problem. 
ies virtual photon is expected to interact with one nucleon 


Smiy. But the interaction amplitude of an electron with a 


bound nucleon is unknown. The natural thing to do is to 
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eeoume that since the deuteron is a loosely bound systen, 
miemOnovon interacts with a nearly free nucleon, while the 
Sener nucleon participates just as a spectator. This is 


called the impulse approximation and is described graphic- 


(4 


Figure 8 


mel by Figure 8. 





Here the bubbles represent the photon - free nucleon vertex. 
Distortion of the free electromagnetic form factor (vertex) 
due to binding, i.e., the difference between bound and free 
BoomeracvOrs, 1S thought to be less than 5% [33]. Now that 
moe impulse approximation has freed the nucleons, the nucleon 
Pemcicipating with the photon in the interaction is described 
by the free field Dirac Equation. 

Since the impulse approximation involves the inter- 
meulon Of a photon with an in aiwslctee "free" nucleon, the 


terenty +15 asain of the form coe. where 0! is siven by 
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u De 2, Uv 
OF = Ff + F O 

3 een apa doa 
WoL Ce Oy _UV 
eae inc yy + Noeds)o  q 


The deuteron current is then 


7 
Ja 





Wonuer Oly. 
p n 


Nee 


= 2 eee) 


C2 ike) 


Adier [45] writes the relevant part of the S-matrix form in 


the impulse approximation as 


2 ta 
o~ [S.,1° ~ IM.,1° = 


Te elt 


ieeceas in the proton case, era i 


to be 


2 
€ 


—_ Eee 


where n = 9°/4Me, and Gy» Crema a 


Q 


[1 + 2(1-n) tan 


M 


(220) 


s evaluated and found 


2 


MIP 


](2-114) 


are the charge, quad- 


rupole and magnetic form factors given by (including their 


meotac Lamits). 








GF (Gy + Ge ) D, a it G22 is) 
Pp n 7 
anda 
2 
= eee a 
Go Ge + iF ) Do MQ (2 aloe) 
p nN 
and 
E M Ma 
Gy = (Gp + Gy )Dyy + (Gy + Gry )2Dyy mw Ug 
p n p n p 
(Qe) 
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Mimere the D's are given in equations (107) through (110). 
ems veps from a point deuteron to a finite size deuteron 
meevhus carried out by the Aer omuction of the free nucleon 
Porm facvors. This step is justified by the impulse 
approximation. 

The Pencributlons of a6 (charge), oR (quadrupole) and 


cr (magnetic) to 


2 e) Gee ree 2 eee 
= + — + = - 
A(Q*) = Goat) + g n°Gglae) + SnGylae) (2-118) 
are shown in Figure 9. 
ui. Neutron Charge Form Factor from e-d Scattering 


Equation (2-117) can be simplified considerably by 


Secimne, 1.€., 


G 


' i= 
WG o » Ug aya 


We d 
and 
“M "M G 
te a as 
p rh © 


Cre ec CG... yu. Dr (2-119) 


Spervumaw With | =~ 7 > 1, 


2 2 


2 12 
[Dg(1 + § muy + 


: 12 7 ee siege 2 : 


= (6. +c, ) ee 2g 9% 


p nN 


Co| = 


ye ple Cae 20) 
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Contributions to A(a*) 


So Quadrupole 


Magnetic —” 





ica. 
Figure 9 
: 2 2 2 2 
The functions oF (charge), 9 n Go (quadrupole) and 3 nGy 
(magnetic). The smooth curve is the sum of these three 


wong Celene iroem Ret, 446). 
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where FS is called the deuteron structure factor and is the 
same as the result of Jankus for a point (nucleons) in the 


@euceron. From (2-120) 


1l+——— |. (entel) 


The aim of the experiment is the evaluation of the 


woweee worm factor of the neutron. The magnetic contribu- 


walom LO ae and Ge must therefore be eliminated. In the case 


Fn 


d 


of the deuteron the magnetic contribution is given by 





120) as 
mene ole bane = (ete) 
5 d 2 
so that 
2 Beye 2 Cox 
ean Coy =P (221,235 
d W@ie te joel 


where the quadrupole term has been neglected. 
In the case of the proton the magnetic contribution is 


found by considering 


ee CG 
qe = eps n _p (1+ 2 bane as (2a) 
p i 1=-n 2 
Scaling assumes Gy = UG, where wu = 2.79 So, 
Gr 
2 2 2 8 
= - + =) ye o-12 
Gy = TS (laren ( 5) 
thie term 
Cee lt? tan’ = (aeeD 


2 


fal 





1s then the magnetic correction term. Thus 


2 
GQ" = 1 + 
ir ( Ce 
Somunat 
2 
G" (1-n) 
a6 5 Sees ey 
p ae s€ 
where 
D O 
Cea ee 
PO nott 


Finally, for experimental comparisons 


"Eg "E 
Dp Dp 
or 
"EL il a 
ae a Re Cane Oa a (2226) 
iB Gieee ok, 
p Dp 


The value of G, from (2-128) is very sensitive to the 
uncertainty in F, and Ge [Gr . For example a 1% error in 
F, translates into Beit a 10% error in Gp /Gp . 

The only reliable data in the range tae oanener 
here have been presented by Drickey and Hand [1]. Their 
results were in complete disagreement with the neutron- 


ZectmOimumveraculom slope, Their data will be presented 


manor ceantalveus Of the theory in part B5. 
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>. Review of Relativistic Deuteron Theories 
Essentially, all the material above is based upon 
the solutions to a nonrelativistic Schroedinger equation for 
the two-nucleon problem. The results produced using the 
nonrelativistic S and D state wave functions avenmal vtames . 
Surprisingly good. Still, there are problems as the nagging 
difficulty in the determination of the percentage D state, 
and the theoretical determination of the deuteron magnetic 
moment, not to mention the conflict that is the G1 oso 
this work, i.e., Gp * 0 while dG, /a(q°) ah 
n n ee 

Hence sufficient reason exists for improving the theory, 
as well as the experiment. 

The Jankus theory on which the cross section (2-106) 
was derived is itself one of the major limitations in the 
development of better answers for the deuteron. Gross 
states that the theory gives no indication of what has been 
Beiveout, Or how to calculate COukec vious, whiten mi cht 
iieiude the effects of of f-the-mass-shell contributions 
arising from the fact that the nucleons are bound and not 
fee. Intermediate states, i.e., pion exchange currents, 
could cause further corrections. Thus to isolate possible 
corrections deoviecclaui vistmemtneory should be developed, 
one that would reduce to Jankus theory at low a Correc- 
tions then could be localized and more carefully defined. 
ietoeis 2 difficult order to iM eanemmordate no entirely 
Savlolacvory solution has been proposed, though many authors 
have tried. The most interesting attempts are briefly dis- 


cussed below. 
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Van [34] has calculated, using deuteron wave functions 
femmved trom a Bevhe-Salpeter equation, relativistic cor- 
Mecerons to the deuteron form factor. These wave functions, 
Woeema are relativistic generalizations of Hulthen wave : 
functions (no hard core) tend to make the functional D of 
Equation (2-106) larger at high momentum transfers than do 


meaeis Navinge a hard core. The essential part of Van's 


results is a comparison between his predictions of Do Ge 
ieven’'s notation) calculated with the (relativistic) 
Hulthen wave functions and the corresponding Do (Fon > Van) 
of the nonrelativistic case. The comparison showed 

POD (g°) > PNR an a range of q® from 0 to 9 F~°, which 


Chi ch 


mesults in Fa Betmarlarger than in the relativistic case. 
Recalling (2-128) quickly confirms the fact that the end 
result is to make Gr [Gy even smaller. The increase in 


n 
ie from oe inovuemenped TO the quadrupole contribution 


2 ey 


(43.5% at qo = 1F  “), and to final state wave function 


Merentz contraction (-1.5% at q* al RUS), 

Adler [45] has obtained what amounts to first order 
relativistic corrections to the nonrelativistic impulse 
PeOrOximervion deuteron form factors. His theory is intrin- 
Peecivemonrelatsavyistic and differs from that of Section 
Miepert 82 only by the use of Pauli spinors (for convenience) 
instead of Dirac spinors, and by retaining powers of 


deuteron moments such that the final cross sections are 


correct to order q°/m. The extra terms now found in the 
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deuteron current are incorporated into additive corrections 


GO Go and Go: specifically 


ADLER Ges 9g 2 ath R 
G =-G,-—5- fa ig eee cs. . NG 
0 


C ic ams la|r G c 
and 
ADLER _ 272 Re rai 
G ee eC 
d 
where 
7 jn(ar/2) 6 f 2 45,ar/2)-43,(ar/2) 
AG, = 24 a We eo Gt rere eer (OT: 
0 a 5 0 la|r 
and where 
G = 1G; eG G, = F a5) eats (q°) 
RS R Re. 2 R il: 2 ; 
p n 
a ~2 2 hoe —5 
ieee = 1.0 F where Go =O! AGG amos x 10 and where 
Gy 0,03, Gh es Te hence the corrections are rather 
Too small to effect Fa ey eDiieso Em toe raght direction. 
Further, these corrections to F, do not alter the value of 


Hy place; is, depends only on P and the value of Qa is 


D? 
changed by a bare 1%. In summary then, it is doubtful that 
these Ceiireevleons are very meaningful in a theory using 
nonrelativistic wave functions. 

Adler glossed over at least one sticky problem in the 
above work. He admits the need for a "fudge factor" in 


resolving the kinematical difficulty in the lab frame of an 


energy transfer (recoil) to only one of the nucleons by 
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the virtual photon. Gourdin [35] has made what amounts to 
the same calculation but in the Breit frame where the 


memiciatlum Cranster 1S entirely spacelike, i.e., q@ = aoe 


a MPs titel bLy Was in the transformation of the 
femeelativistic deuteron wave functions into the Breit 
frame. His solution to this was to assume that the wave 
function were unchanged. Gourdin's results were approxi- 
mately the inverse of Adler's, i.e., of the same magnitude 
Dut Opposite in sign. Unlike Adler, Gourdin did derive a 
Tao Correction to the magnetic form factor but it is 


so small that it is not observed for q° 0 ae 


Tne wy ale 
Umemeot bese works lies in the fact that it almost certainly 
seals the fate of any straightforward attempt to determine 
(quasi) relativistic corrections. 

The strictly relativistic approaches (a covariant 
dispersion relation representation) are more basic and more 
adifficult than the quasi-relativistic work above. success 
WiGheuimis vechinigque is limited to the low momentum transfer 
region where it is now believed that the theory is more 
accurate than experiment. The early papers of Jones [4] and 
Cutkosky [36] are interesting when applied to the Gr 
problem, Using the triangle impulse approximation seeean 
(page 67) Cutkosky calculated the simplest absorptive 


amplitude which yields a deuteron structure factor of (all 


spins ignored) 


Yn cot = 
cn a, |= = uM. (2-129) 
. Yn cot M Valen 
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iMemrelabivistic analog using the tail approximation is 


oO 





ae ; 
Fifa") - [ wri gtar/2)ar = 4 ean Jal (2-130) 
A a 
where 
_ Siemens 
T= Mott ME Pas com 
: 7 ae ee 
Comparison of the two Bs and Fa wie idwthe relavivist uc 
Comrect ions 
Fl a EBe) 
d et “13 
x 
-> V 
and |g | ——- — (2~133) 
YV1l+n 


These imply that a surprisingly rapid curvature is needed at 


small ae for Gy in order that the n-e interaction be in 


n 
consonance with G —OmayenO, 3 = 0, © Boe, 


E 
n 


Seoss Nasedone the most extensive work in the area of 
the purely relativistic approach [2, 8, 9, 10, 11]. He 
ambitiously set out to construct a complete theory fully 
Poareeevistae and fundamental in the sense that very few 
phenomenological parameters occur and based on single- 
Vvomlaple Unsubtracted dispersion relations and coupled 
unitarity equations. He had partial success but numerical 
PaoGmerlonmcn Ol sume rau10 Of relativistic to nonrelativistic 
O 


Peounmraceors scems much too large at high q in the: range 


of low q? Cones = 1.0 F7e) agreement is reasonable, and 


his predictions for Ha and Qa aeommiltham 2% of the 
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aaeerimenval Valve. in work on relativistic corrections to 
the impulse approximation in elastic e-d scattering Gross 
realized that the purely dispersion technique tended to 

teem cthesinternal structure of the interacting particles, 
i1,e€., the deuteron wave function and the nucleon form fac- 
tors with the deuteron current. This leads him to eventually 
~emeeos Nis theory in a formulation closer to that of 
potential theory. 

Specifically he obtained a modified Bethe-Salpeter type 
wave function retaining terms up to the order of m7e and 
then used the relativistic nucleon current to derive the 
fameeron form factor. An estimate of the size of the cor- 
Gemwron LO the charge form factor is obtained by taking the 
Buome Of the correction term at q@ “Oe aross found that 
the correction is proportional to ~1/3M* Mine 1S) ef 
Paice sign SO as to Make the neutron charge form factor 
more positive (tcwards agreement with thermal n-e slope 
data). 

The relativistic correction mentioned above came mainly 
Paotimerelativistie Modifications of the deuteron wave func- 
brens and the nucleon current.. Adler had previously looked 
into the nucleon current modifications but his analysis had 
been based on nonrelativistic deuteron wave functions. 
Gress States that the wave function modifications have the 
fapeest ettect and hence Adler's ee oer aons are insufficient. 


Van did use relativistic wave functions from a Bethe-Salpeter 


Vineness TWdyeOoL Corrections to the impulse 
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aeemextmavion. His Gifficulty was in using a separable 
potential which was equivalent to using a Hulthen wave 
function which has no hard core, and thus gives values of 
Ba E@@mlaree., he trouble 1s that here relativistic cor- 
rections to e-d scattering and the model dependence 
question have been intermixed. Van did M@ Peron Osse eb has 9D y 
weemreaon Ol his relativistic corrections to corresponding 
Poreeiwavivisctic models since the two were not in the same 
reference frame. That is to say the relativistic model (in 
the Breit frame for example) when Lorentz transformed into 
PeemceuLeron rest frame was not compared to a corresponding 
nonrelativistic model, which if it is a good model already, 
Mieco lt melaravistac effects as seen in the rest 
eine 

icmp iiicipal COnUribuLcion bo the correction of the 
charge form factor of the deuteron comes from taking into 
MaoommumuiesrecOll MOtiom Of the deuteron, which in turn 
Seoveryus tie wave function. Lorentz contraction also adds 
Pen emave 1UNCb1OnN distortion, though not as much as in 
Vermeewcace, Lenoring the correction for the quadrupole and 
lapowe@rc terms, the result of Gross' calculation for the 


charce form factor is 


Go = (Gp + G, ) 
p nN 


R 
Do (2-134) 


*Van did this in order to improve on Jones and Cutkosky's 
ear 





asymptotic "tail" wave function approximation ~, = => 
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where 


- n Bye 
ae = aad f (uo +w') J, aie AD. (2-135) 
0 


and 


ai a 1 a ae , ‘ 
AD«, af uru uku + 7w wrw wkw) (jot55) dr 
0 


+ nf (uru' + wrw') Jj, ar. (221346) 
0 


In this expression prime refers to differentiation with 
iteoeecy FOr, and 


KE re a ~ a*l. (Bis 
aie 


The arguments of the Bessel functions are of course qr/2, 


emeec 1S calle emi wets Ol sGrOoss 'S COrrections , A(q‘) is 


A(Q°) = (G. +G, De = (a, +G, )<(FS + AF ,) (2-138) 
p n p n 


where the quadrupole contribution to the charge form factor, 
8/9n°G® , is small. Thus, since A is an experimental quan- 


Pi vercoma dees Nou Change due to corrections, a change is 


iereed in Gy » To see this write 
n 
us NR,.NR _. NR = 
Go = (Gp + Gy )Do - (Gp + Gp )Do (2-139) 
p n p nN 


(no corrections). 
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es a NR 
Go = (Gy, ae )D, = (G te [G, 


+ AG, 1) (Dan 7 ADQ) (2-140) 
p n p n n 


CwMterelaAtivistile corrections). 


Equating the right hand sides of (2-139) and (2-140) yields 


NR . NR 
-AG, Do = (Gp + Gp) AD, =~ G, AD, 
n p n p 
Or 
ADo AF 
Cc (25058) 
E, E, pur ee 


where Bg Pee macuLbeton’s Charge structure function. Gross 


wrote (2-141) leaving out the F Since F, < 1 for very low 


q° a 

Gis AG, ~ Gp AF, Sives an Upper) limit to the correction as 

Cee 0 ie ae interaction. Figure 10 summarizes the 

femrecugens, in the range of qgé Con smeemead here. Gross" 

correction AG, homcOopory tonal, to — 25 anc Bagot ica. ly 

model Reo snaent.. = ct a 
womioeenis —eOrrection to the charge form factor of the 


Memeroneplus Feshbach—Lomon wave functions for the deuteron, 


Gross reanalyzed the Gr =“ 0 data of Drickey and Hand. 
n 


Figure 11 shows a comparison of results with Partovi (i.e. 
HJ) and Feshbach-Lomon wave functions. Gross! correction 
coupled with FL wave functions seems to remove the discrep- 
Bae mDelweemermne electron scattering and Hee Mireracy Lon 


results. 
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/ 
y, D AGp-(Partovi) 
0.03 ore q 
/ eae 
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/ oe 
0.02 ys = 
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yi ee 
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OU | i ee 
Tg aGe (Partovi ) 
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// ee 
iy —————) 
0 eG) a 10) Be) Hg pea) 6.0 fen) 
Figure 10 
The relativistic correction is labeled by asa the 


En 
difference between a Fe Bate Lomon model and Ghee wrarcovs 
model is labeled Gprh ah) and the meson exchange contri- 
butions are labeled. sok : “The theemaleneuuvron Slope is 
given for reference (eaen feom Ret. 2). 
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ae) 





PEieure 1 1. 


imemdabvalet Driekey and Hand. The dashed lines are the 
thermal neutron slope, +0.0193 + 0.0004 and the fit of 
Chilton and Urhane (unpublished). The number on the right 
Demo wc treure 25 the slope determined from the first three 
data points of each graph. (a) The data analyzed with 
Partovi wave function, (b) the data analyzed with Partovi 
Ti mointrronmimcindins relativistic corrections, and 

(c) the data analyzed with the Feshbach-Lomon wave function 
including relativistic corrections (taken from Ref. 2). 
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Pere cetiit auucmpoul bY Chailton and Urhane [37] to explain 
aemee 0 for q° = 0,3 ~ 1.0 PF“ 
n n 
Daeroeine che nucleon form factors with rational functions 


with Gy (0) # 0, was made 


of the variables 1 + a°/m,*. They found Gy, is never far 
i 
from zero and bends over rapidly with increasing qt (see 


Figure 10). Nevertheless with Chilton and Urhane's fit, 


G, does not get back to zero till beyond 12 R7e, 
n 
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een SAND PERIMENTAL PROCEDURES 


Peel NPS LINAC 

The Naval Postgraduate School linear accelerator (LINAC) 
is a low intensity (20 uamp maximum average current), 100 MeV, 
three section, 30 foot long electron linear accelerator. 
While its construction and general operation have been 
described elsewhere [38], several important modifications 
have been accomplished and will be mentioned briefly. 

A supplementary cooling water system with a steam heat 
exchange was added for the accelerating sections. This was 
Oe so that the water temperature could be precisely con- 
trolled (41°C) throughout a run lasting for an extended 
length of time. In the present experiment, runs of up to 
38 hours were made, and since the machine's energy stability 
meee oconely dependent on the accelerator's changes in 
Piieetravure, Chis 1S a critical point. 

Uieeeceny ly tne Garset chamber vacuum had been 
isolated from the accelerator sections by a three mil 
aluminum window 32 inches upstream from the target. This 
Window caused considerable beam spread and by the time the 
target was illuminated by electrons, the beam spot was 
approximately 3/16 - 1/4 inch in diameter. Now, with the 
addition of a fore-diffusion vacuum pump combination in the 
Careet area, the entire accelerator, from electron gun 


See@urecy sto cerfet. 25 Open. In fact with the spectrometer 
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PooveastO une target chamber, electrons pass only through 
miemoarses itself and a five mil aluminum exit window from 
Che spectrometer before striking the detectors. Upstream 
secondary emission monitors (SEM) have been moved to loca- 
tions downstream of the target. The accelerator - target 
chamber - spectrometer combination etn eainied at a 
pressure of 107° mm of He. The removal of the SEM and 
Peeagow nas Significantly reduced the beam spot size, and 
it is now estimated from beam appearance on a zinc-sulfide 
(ZnS) screen that 2/3 of the electrons are confined to a 
spot 1.0 mm in diameter, 

Past nuclear structure experiments [39] have experienced 
Severe problems in background radiation levels. Major 
shielding modifications were made to correct this situation. 
More paraffin and borated paraffin in the form of a floor 
tO ceiling wall was added to the existing shielding in the 
beam deflection area. In the target area, a platform loaded 
fee aratfin four feet above the floor and extending 30 
inches out from the wall was suspended from the overhead. 
Its purpose was to shield the counter house from line-of- 
sight radiation from the beam dump. Wax was placed in the 
spaces between the IJ-beams in the overhead in order to 
reduce reflection of neutrons into the counting system. 
iimoriiawywslitelding effort was on the counting house 
itself. Concrete blocks were replaced by boxes of paraffin 
and much lead was replaced by special bricks composed of 


Bb bere leaceby welpht end polyethylene. <A cross 
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Secvion of the counter house is shown in Figure 12. This 
shielding arrangement reduced the background to quite an 


Meeenvabvle ratio of 2000 to l. 


B. THE TEN-CHANNEL COUNTING SYSTEM 

fae oc leon—-Channel System 

In general there are two main categories of detec- 

tors: slit-defined momentum acceptance, and counter- 
Pemm@etry defined acceptance detectors. A single counter, 
slit defined detector was replaced by a ten counter 
Pewermry detined system. The ten counters, which span a 
portion of the focal plane of the spectrometer and define 
adjacent momentum intervals, are known as a "ladder." The 
Pareatares of this arrangement over the single channel 
owes ObVi0US; for one spectrometer momentum setting, 
ten pieces of information instead of one are gathered. 

iiewladder is composed to ten 7/16 inch high, 1/16 inch 
thick NE102 plastic scintillators placed approximately 1/2 
inch behind the focal plane of the spectrometer as shown 
in Figure oe A large one-piece scintillator used in 
Eomm@emdence wich the front counters and called the backing 
PeuNverwsic placed 3 inches behind the 10 singles counters. 
The choice of 7/16 inch high (in the dispersion direction) 


counters produce a momentum acceptance Ap/p of about 0.3% 


*The ten scintillators (counters or channels) are num- 
pecimiemeaseendineg order 1 to 10, where channel 1 represents 
the lowest momentum accepted. 
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per channel or about 3% for the entire ladder. An electron 
aeamectea from the target and subject to the spectrometer's 


meaemervtce tteld B, follows a path of radius r, given by 


ee eee x 
<6 Bec — (3-1) 


Thus higher momentum electrons are seen by the higher 
numbered channels. Further the momentum resolution is 
greatest for the higher channels. An electron traversing 
Weiemopecurometer and passing through a particular scintil- 
lator produces a pulse that is collected by an Amperex 
Paeeo pnovomultiplier tube, amplified, transmitted, 
Pereiamiaaved, fanned out, and finally registered on one 
Gmmmenmescalers. A block diagram of the counting system is 
shown in Figure sie With the exception of the scalers, 
Wiomemeallow a2 counting speed of 20 MHz, the electronics is 


capable of 100 MHz counting rates. 


( 


EeeococOlur on and Eificiency of the Ten-Channel System 


In order:-that the system be useful, two quantities 
must be known: The relative counting efficiency of each 
channel and the momentum separation of adjacent channels. 
The momentum separation (resolution) of the channels was 
Savcened by taking elastic scattering peaks from a at 


carbon foil and detecting the electrons scattered as a 


*The Carovevpeecowana BG and G counting Equipment was 
loaned to NPS by the High Energy Physics Laboratory, Stanford 
Wnavers iy Vy. 
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function of spectrometer momentum in each channel. That 
is, the 10 counters were treated as if they were single 
counters while the sharp elastic peak was moved through 
them by changing the momentum setting of the spectrometer. 
A plot of counts in the scalers vs. spectrometer momentum 
yields 10 elastic scattering peaks (Figure 14). The separ- 
ation of the peaks as a function of electron momentum p is 
a measure of the channel's resolution. The 10 elastic peaks 
Muetecd are not identical in height or area, but since 
they all "see" the same elastic carbon peak, the differences 
are ascribed to the differing relative efficiencies in detec 
tion and recording of the electrons. 

Mowespecifically, the resolution of the 10 channels was 
determmined by taking 21 - ten peak sets at incident elec- 
tron momentums from 27 a to 94 “eV and with scattering 


angles varying between 45° and 105°. These experimental 


Gemeewere Tit CO a quadratic expression 
a Ze 
p = p.(a + BN + YN‘) (3-2) 


where N is‘’the channel number, Pe is the momentum setting of 
the spectrometer, and p is the mean momentum of ;whanne lL N. 

Pee eoerileeprodueca values for a, 6 and Y. If oe is defined 
as the momentum setting for which a particle of momentum P, 
entering the spectrometer passes through the center of 
channel 6, then (3-2) can be rewritten as 


s 7 ye 
Dy (P,) = eet ea 6) +B(N-6)"] 
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A and B again are constants independent of p. The channel 
momentum acceptance is then 


d(py(p.)) 
ainasy/ Ps = A + 2B(N-6). ea 


itvrcal values of A and B are 3.0 x 107° ana ein 
Mesepecctively and it 1s clearly seen that the higher channels 
jem N > 5) have a higher resolution than the lower 
@m@aanels. 

So that the energy spectrum of the scattered electrons 
can be determined properly the relative efficiencies need 
to be measured. The efficiencies can vary for several 
reasons, some of which are: Differences in detection 
bibesmolds and physical properties of the plastic scintil- 
Powers wopecuromecver transmission characteristics; and 
solid angle. The determination of the relative efficiencies 
was based on the method of Crannel and Suelzle [40], and 
essentially involved measuring the shape of an inelastic 
continuum of scattered electrons. The shape was found by 
observing the spectrum with all ten channels at several 
different momentum settings assuming initially that each 
Penne has efficiency one. Then the counts from individual 
channels are compared to the calculated shape and corrected 
channel efficiencies are obtained. Next, a more accurate 
ep mr omorecameved on the basis of the corrected ea tae mencies 
and the process is repeated. Three or four iterations will 


normally produce constant efficiencies. 


oS 





ivparecalily, tive different momentum settings over a 
total momentum range of 6% (twice the 10 channel span) with 
at least 1000 counts per channel per setting were taken. | 
It is important that the efficiency data be taken in the 
energy range in which the experiment 1s Como Der Ormod. 
meme@enmtt an clastic deuteron peak is located at 50 MeV, then 
the machine energy should be adjusted so that a smooth part 
of some inelastic continuum lies at that point. In carbon 
ieee radiation tail of the elastic peak is sufficiently 
smooth at 30 MeV below the peak, thus for the above situa- 
tion a machine energy of 80 MeV would be set, and the five 
momentum (energy) settings could be taken at 52, 51, 50, 
49, and 48 MeV for example. Carbon is obviously not suited 
for the NPS machine because then @ practical limit of 
65 MeV would be imposed on the incident energy for the 
proton and deuteron peaks (i.e., 65 + 30 MeV = 95 MeV 
needed for efficiency and any more would exceed machine 
Poepobitsicvy)., Aluminum was substituted; its inelastic tail 
Pompoincwpeaks Only at 2.2 and 3.0 MeV. Arbitrarily a value 
of.6 MeV down from the aluminum-elastic peak was chosen as 
the experimental area from which to obtain the relative 
counting efficiencies. Normally it is best to exactly dup- 
ittceereminebne Ciliicienecy measurement the conditions of the 
particular experiment, but this was found to be impossible 
in some cases. On the high energy runs, 94 MeV (location of 
elastic peaks under consideration), 94 + 6 = 100 MeV was an 


impractical energy to obtain and maintain. Thus the 
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efficiency data were taken at a somewhat lower energy. Also 
aluminum not CH, oie CD, was always used for efficiencies 
pees t Was desirable not to take the risk of burning or 
Peomroying the experimental targets during the long high 
Peamesurength used in efficiency runs. 

To calculate the relative efficiency of channel k it 
1s necessary to predict the total number of counts channel 
feeeecives in the L spectrometer momentum settings. Thus 

LF L Py, (KP, ) 
ete.) - > E.(p,p.)T(p.)dp (3-5) 
s=] s=1] p,(k,p.) 
where p is the particle momentun, T(p.) dp is the number of 
Sweemeols With momentum between p and p + dp which enter 
the spectrometer, E,(p sp.) is the efficiency of channel k 
for particles of momentum p, and Py, (kP,) ana Dig laabee 
the highest and lowest momentum accepted by channel k. 


Assume that 


2 

- alate: 

T(p,) = ) RP (3-6) 
n=0 


where a can be determined (needs to be done only once) by 
10 Ti 

a least squares fit to the points ae = mee C,(p.) as 

aeeuineu.On Of De: Hence the efficiency of channel k is 


prvermupy the ratio 


L L E 
deeds 
E.=| 5 ¢,(p.)/ > > Raps (3-7) 
s=1 s=l1 n=0 
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which is the total number of counts observed vs. the pre- 


dicted total number of counts. 


See LHe TARGETS 

Solid polyethylene (CH,),, and deuterized polyethylene 
(CD,),, were chosen as the sources of protons and deuterons. 
igemrmmediate difficulty is that these plastic targets would 
be destroyed by heat (melted) in the concentrated electron 
beam. Hence the CH, and CD, targets were made to move in 
the horizontal direction so that the beam sampled a width 
of about one inch over the target face. The carbon in CH, 
and CD, has much greater atomic weight (12 umu) than does 
Emeier 1s5otope of hydrogen (H = 1.00783 umu and D=2.01410 
umu) and hence suffers much less recoil when struck by an 
electron. Thus both the elastic proton and deuteron peaks 
ride well down on the inelastic tail of their respective 
carbon peaks. In order to isolate the proton and deuteron 
peaks it 1S necessary to subtract out the carbon content. 
A carbon target is needed then for standardization. Hence 
arrayed down the target ladder were movable CD, and CH, 
targets, a stationary carbon standardization target, alum- 
inum target for efficiencies and a ZnS screen for machine 
Pian and beam observation. 

The CD. target was 3 cm x 10 cm x 0.10122 cm and 
specified to be 99.9% enriched. The thickness of the CH. 


and C targets were governed by the desire to make the number 


Of C atoms in each of the three targets as nearly the same 
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as possible. With the material available a CH, target 
ieorst + 0.00070 cm thick and a 0.09398 + 0.00061 em thick 
Carbon target were constructed. The carbon target was 
produced by laminating 7 thin (5 mils) graphite foils. The 
densities of the samples were determined to be: p(CD.) = 
moor +t ,006 gm em >, p (CH) ero ese 02005 om em, and 
o(C) = 1.028 + 0.005 gm em™>. 

In the course of a run the electron beam illuminated 
the CH, and CD. targets for as long as 8 hours apiece. 
After these long runs considerable browning of the material 
Pema rocnic, and the possibility of damage to the target 
because of a change in the composition of the material was 
considered. Subsequent runs on the undarkened area versus 
bie Garkened area gave no indication of any decomposition. 
Piieemer the darkened areas tended to fade in time so no 
permanent damage could be supposed. The browning of the 
target iS assumed to be an atomic electron dislocation 


phenomena, 


D. COLLECTION OF DATA 


Pepe rececioing Upon the momentum transfer to investigate, 


the appropriate incident energy-scattering angle combination 
was chosen. The basic measurement desired was the ratio of 

the electron-proton to the electron-deutron scattering 

eross section. Forward scattering angles were chosen where 

EOssmole In Order Lo emphasize the electric interaction. 


Withecmiis eaneular restriction the momentum transfers ranged 


ois 


4 
am 





2 @ 


emete provon and deuteron from about 0.1 < q~ < 0.5 F- 
for the available machine energies. 

The machine was tuned for maximum beam and current 
Seeotlaty. The energy defining slits were set for 
AE/E ~ 0.5%, and the beam spot was carefully steered to the 
mieeet Center line, Collection of raw data consists of 
recording the number of counts shown in scalers 1 through 
10, the backing scaler counts, the time elapsed while 
counting, the spectrometer setting in MeV, and the number 
Siecilectrons passing through the target. All the scaler 
Gaver are fed directly to a teletype machine that both punches 
an 8 track paper tape, and prints out a hard copy. The 
paper tape is then interpreted and IBM cards produced for 
use in computer data analysis. 

For this experiment data were taken in the following 
way. The machine was tuned to an energy approximately 6 MeV 
higher than the energy at which the elastic proton peak 
mommecmoccin, Eiiiciency Gata using a 20 mil aluminum target 
were taken. The machine energy was then lowered to the 
value selected for the proton run. Comparison elastic 


carbon peaks in CH, and C were taken to provide normaliza- 


2 
Mihewmeaba Used in Scaling the carbon for subtraction. Then 
the elastic proton peak from CH, and the underlying carbon 
tail for subtraction were obtained. Machine energy was 
again lowered to the proper value and the deuteron runs were 


eomplevecsam a Manner Similar to that for the proton. On 


Seleawse the COUnving rates in the front channels were 
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immived to (5 counts per second, while the backing channel 
was limited to 100 counts per sec. Tests showed that the 
counting system (presumably the scalers themselves, not 

the counting electronics) could digest counting rates of up 
to 200 counts per second before any appreciable loss was 
Semwected, Since the deuteron is eee as heavy as the 


proton, less incident energy is required to achieve the same 


momentum transfer as that of the proton. In all experiments *& 


a peak was determined by three spectrometer settings, hence 
30 points define a scattering peak. The small energy steps 
taken depended on the machine energy and the relative energy 
(momentum) separation of the various channels. The aim was 
Pemachieve a uniform distribution of points over the area 
investigated. 

Hiiemomunespedk 15 located on the ladder is important. 
iMitiemena channels, 1, 2, 9, and 10 tend to have more widely 
Melee eetiveiencies than the center channels and they 
Seem temee Cage effects, so they are unsuitable for the 
peak determination. It was found that placing the peak 
initially between channels 6 and 7 gave the best results. 
Then the top of the peak was defined by a single counter 
reducing a potential problem of matching the counts of 
ft terenmechannels av a critical point. The peak height 
Was) chosen on the basis of time versus statistical accuracy. 
itcvweres. a total of 4 x 10" counts are needed under a peak 
to insure 0.5% counting statistics, but the time needed to 


meeMmilmeate that Many counvtvs 1S Often prohibitive, 


oe 











particularly at the higher values of ae It was decided 


meateno more than three hours per spectrometer setting 
Pemrasbpe allotted for counting. 

itmeerder to extend the range of momentum transfers, 
another set of data was taken by the author at the Mark III 
accelerator at Stanford University. Where possible, over- 
lapping sets of data were taken to check for any systematic 
errors in the equipment. The range covered at Stanford 


£ e No systematic deviations between 


meemo.20 < q < 0.80 F 
Pimewo sels Of measurements were found and the data 
Power, ol was performed independent of the source. 

The experimental equipment at Stanford is basically the 
same as that at the Naval Postgraduate School with the 
exception that a 100 channel ladder is used and that part 


of the data reduction is performed on an IBM 7700 on-line 


Somourer, 
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i. SOA RE DUCT TON 


A. ANALYSIS OF RAW DATA 

havecdava are accumulated in the form of counts in ten 
channels for a single spectrometer setting, and are analyzed 
mame olded' by matching the counts in a particular channel 
to the appropriate energy. The assignment of energy values 
to channels is determined by (3-3). The raw counts are 
first corrected for background, then counting rate correc- 
tions are applied. Finally the efficiency determined by 
(3-7) is used to predict a corrected count for the partic- 
Wiereecnannel, In practice all the above calculations are 
Carried out by computer. Typical unfolded energy spectra 
fOr both proton and deuteron runs are shown in Figures 15, 
ome y, and 18. 
B. CALCULATION OF EXPERIMENTAL CROSS SECTIONS AND FORM 

FACTORS 

iiemwexperimental cross sections for e-p and e-d scatter- 


ing were determined from 


p,d N 


io sc Pea -] a 
rae = Ny, n, ao KK, em” ster CAs) 


where ne, is the number of scattered electrons detected by 
Ciemwectromever, NN is the number of incident electrons, 
in 


feereene number of target atoms per om”, ieeke, bie Solad 


t 


angle subtended by the spectrometer entrance port, and 
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finally K, and K, are the radiative corrections of Tsai rae) 
and Bethe-Ashkin [42]. 

At the heart of the calculation of (4-1) is N., which is 
the integrated arvea under the scattering peak in counts-—MeV 
per MeV. The scattering peak is produced by drawing a smooth 
eeevesunrough the data points in the counts vs. energy plot. 

Both the elastic proton and deuteron peaks ride on the 
radiation tail of the elastic carbon peak (see Figure 15, 
ete.). Thus, in order to obtain the relevant areas for 
the cross section, the underlying carbon tail must be sub- 
tracted out. Comparison of the elastic carbon peaks in 
CH, and CD. to that of the pure carbon target provided the 
target normalization factors. These factors consisted of 
Meme chatts’ of the elastic carbon peak from the carbon 
target such that its peak energy coincided with that of the 
carbon peak from CH, or CD,3 and (2) vertical "scaling" of 
the carbon tail through a comparison of the areas of the 
elastic carbon peaks. The area remaining after the sub- 
meeeo 4 ON 1S Nee 

The number of electrons incident upon the target is 
measured by a Secondary Emission Monitor (SEM) downstream 
from the target. The SEM is not an absolute device to 
measure the charge. In a ratio experiment, the only 
requirement is that it be stable during a set of runs. 
This stability has been frequently checked against a 


Faraday cup. Obviously the number measured is wrong by at 
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ieee the number of electrons scattered by the target, but 
ps 1S SO Small compared to Nae eemco oe LOLally insig-— 
feteLcant. ny is a static parameter; it depends only on 
target composition and target angle $6, 
ot N 
a 
———— , ae 
t A cos $¢ ( 


3 


where p is the target density in gm cm ~, t thickness in em, 
No Avogadro's number, A the atomic weight, and @o the angle 
between the beam line and the downstream target normal. The 


solid angle AXis defined by 1 inch thick tungsten slits on 


the spectrometer entrance port, 16 inches from target center. 


It can be seen by comparison of two experimental cross 
Peertons that the only remaining quantities that enter a 
ratio experiment are Nae and nhs 

As mentioned above, we have used an expression due to 
meeeeoraccouny for the radiative correction Ke: aS C Ola 
rection is a significant improvement over that given by 
Schwinger [43]. Presented here is Schwinger's much simpler 


expression (4-3) as an illustration of the parameters 


Times thie correction. 


a 
6 ORM 
Pa E les Cc il ey 
fee, 6, = )}!en —— - 2 ne O-2)+ 2 
S S Tr | AE £3/° || — Fo eS 36 


(OU (4-4) 


i cos  &n2 AE Bole 


LOR 


{ 





ipieaefiy , Ne iS a multiplicative correction to the cross 
section which accounts for radiation emitted during a large- 
angle scattering event (nuclear pronest rane) . Ky corrects 
Rot ene radiation that occurs before or after the large-angle 
event and is sometimes known as thick target bremsstrahlung. 
For a discussion of the terms in (4-3) and (4-4) see Gordon 
[44]. The § used above corresponds to the nucleon recoil 
factor in (2-77). 

The ratio of the experimental cross section to the cross 
section predicted for a point nucleon (Mott cross section) 
Equation (2-76), yields the experimental form factors ae 
Or ae Me these quantities, the magnetic corrections as 
outlined by equations (2-122) and (2-126) are then applied. 

This yields the square of the electric form factor of 


Pile protlon , ee mon thes SUNVOLetie Squares of the electric 


E 
P 2 8 #2 42 
and quadrupole form factors of the deuteron, Ge 4p 9 n Go : 
d d 
Because of the smallness of n in our range of momentum 


o = C250 a the quadru- 


y 


pole term contributes only about one part in 10 ° to the 


- ae 
transfers Cae peewe x 10 at gq 


form factor and can be safely neglected. 


The end result of the experiment is then the ratio 


Ge /Ge from which the ratio Gp /Gr is extracted with the 


d p n p 
help of equation (2-128). 


In summary, the quantity evaluated from the data is 


2 

G d p d p d d 

eee te Mott * Ss * Sb (4-5) 
Pei d p d p p 

“E,, ea ye ny xe. Xx Ovott x Ke X Ky 
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eee © represents the magnetic correction, the Mott 


“Mott 
Cress section, and the indices p,d refer to proton or 


deuteron, respectively. 


C. EXPERIMENTAL ERRORS 

Brrors encountered in the measurements are of both 
Meeetecvical and systematic nature. Systematic errors, such 
wmedemvalte Of the solid angle, counting rate correction, 
target angle, scattering angle, and SEM efficiency cancelled 
out since they were measured under the same experimental 
eonditions. The absolute error in the energy of the incom- 
pmeeelectrons is about 0.2%. Since the Mott cross section 
for the deuteron and the proton are evaluated at nearly the 
Same energy, this error becomes insignificant. The same is 
true for the errors in the radiative corrections if they are 
evaluated over roughly the same energy cut off AE (Equation 
ao) 

The only significant error of a systematic nature is 
Meemerror in the target thickness and density, pt. Because 
Periewriianness Of the targets, it could be established to 
Omily 0.5%. 

ipemabapisipical error consists of two parts, the actual 
fodiens error and the error in the efficiency correction. 
ims an principle possible to keep these errors arbitrarily 
Piel by taking a very high number of counts. Restrictions 
in counting speed as posed by the electronic equipment and 


Pi MmMMmcrat Lom in the wee of the accelerator pose a 
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practical limit. All but two peaks contain 40,000 counts, 
which would make their statistical error to be 0.5%. The 
Be ouered carbon subtraction introduces an additional en 
depending on the number of carbon counts under the peak area. 
ii eoweor aievne efficiency corrections is quite ini- 
formly 0.5% for the Naval Postgraduate School runs, where 
Povuephily the same number of counts have been taken to 
establish the relative efficiencies. In the Stanford runs 
the error varies from 0.27% to 0.65%. 


For each measurement of ae /G6 : 


ae 
eameresroct, Of the quadratic sum of the counting, efficiency, 


the error is then the 


paommGareet error. 
One half of this percentage error enters the ratio 
Ge. [Gy poem unererore, the calculation of Gy [Gy Sees 


ad =p ae 
is the error quoted in Table VII. 
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V. COMPARISON TO THEORY AND DISCUSSION 


Pee ALE RIMENTAL RESULTS 

The experimental data are presented in Table VI. The 
eamoreawoved consists of counting errors (including the 
carbon subtractions and efficiency corrections) and an error 
Sumvroe in the target thickness. 

Table VII shows four sets of values of Gp [Gr as 
obtained by: (a) the Feshbach-Lomon wave Ree with 
relativistic corrections; (b) the Feshbach-Lomon wave func- 
tion without relativistic corrections; (c) the Partovi 
wave function with relativistic corrections, and; (4d) the 


Partovi wave function without relativistic corrections. 


feeuean Or the relativistic correction AG,, and, as an 
n 


Soil w ene Structure factor fF. for the Feshbach—-Lomon 


d 
wave function are also given. In each of the columms are 
shown the slopes (dG, [Gy ) /dq® to illustrate the rapid 
decrease in slope nee ee Pronlmarelapavists1cal ly 
corrected Feshbach-Lomon wave function to an uncorrected 
Fagpovi wave function. 

iim nnit @f de Vries [28] for the absolute proton 
charge form factor was used to extract Gy from the Gh fe 

n n 


ratio. This fit is in very good agreement with the absolute 


measurements of G by Drickey and Hand [1]. Table VIII 


E 
p 
wae mie ereslilts ot phe b' fit for the two extreme cases 


presented in Table VII. The slopes as determined by the 


dligtsclt 





the two sets of data are also given for comparison with 


the neutron-electron interaction slope. 


Pee YP LSCUSSION 

Mine Figures 9 and 20 illustrate the results tabulated 
in Tables VII and VIII respectively. It can be seen from 
the slope of Figure 20 that the Partovi wave function with- 
ememrelavivistic corrections does not significantly differ 
from the data of Drickey and Hand pictured in Figure 10(a). 
iat is, they bo®h essentially agree that the slope is very 
small or zero. On the other hand the slope for the Feshbach- 
Homer wave function with relativistic corrections approaches 
closely the value of the slope for the neutron-electron 
Mgeract lon. 

Barring unexpected fluctuations in the proton form 


Hercmoer G in the range of momentum transfers of interest 


E, 
D 


here it can be concluded that: 

(a) The Feshbach-Lomon wave function together with the 
Zao vistic conrections AG, MEGCIIOVMeS ste wapparens Giscrep— 
ancy between the eee Slope au q® ~ 0 and the 
slope given from values of oe (obtained by electron scat- 
tering) in the range 0.10 < q@ < 0.80 cae 

(b) Even within the relatively large errors propagated 
aA © Gp , the Partovi wave function with relativistic 

n 


corrections applied, is in disagreement with the neutron- 


plecueemeinteracvion slope. 


I aa 





(c) With the present small experimental uncertainties, 
ie is reasonable to state that this experiment should be 
added to the group of experiments which must be explained 
meen aeuceron theory. That is, a proper model for the 
Geuteron must be found that will predict the slope as 
devermined by the neutron-electron scattering experiments. 

The important qualitative conclusion from these experi- 
mental results is that the neutron has a non-zero charge 
hemmmeiacuor, This implies that there is a charge distribu- 
BLOM Wictvhin the neutron. The neutron thus must have a 
ievetm@ea charge structure with the outermost part of the 


charge distribution being negative in sign. 
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TABLE Viit 


pLOPE OF THE CHARGE FORM FACTOR 






2 
dG. /aq 
n 
2 
q p Grv+AG, 
(F-©) | de Vries nN n 









etal 








0.9881 -0.0027 | -0.0020+ 0.0054 






0.9765 iO neu 30 





OF O0 si 





0.9652 OOOO m at One One meee 10) 15) 






0.9540 tO. 00 7 OCs 






+0,0067+ 0.0074 





0.9431 eOrOany 












Ooi sya al #0.0011 | 40.0047 = 0.0066 
















+0.0184 + 0.0059 


cee eee ee 





Oro 115 +0.0138 





i ee 


-0.0038 + 
Teen) 17 = 
-0.0037 + 
+716) (OMG) sie) = 
-0.0015+H 
-0.00H42 + 


+0.0076 + 


oe 


OF 


.0054 
U0 5 
V00s 5 
PO 3c 
.0074 


0066 
0059 


0-0026 = 0.0036 


Slope of n-e interaction = 0.0193+ 0.0004 


a: No error has been applied to Gp 
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APPENDIX A 
THE NEUTRON CHARGE FORM FACTOR FROM THE 
NEUTRON-ELECTRON INTERACTION 
In the n-e experiments what is measured is a scattering 

amplitude of a neutron by a bound electron. In the Born 
approximation for neutrons of very long wavelength this is 
related to the volume integral of V(r), the interaction 
potential [25]. The magnitude of the n-e interaction is 


expressed in the form of an equivalent potential V, over a 


0 
2 
sphere of radius rg = = oe the classical electron radius, 
m_c 


so that the volume integral of V 





0 is the experimental value 


of the potential. The potential V, is simply a convention 


0 


Since the classical electron radius plays no fundamental 


fever in the problem. Thus, 





O 
free = vo f ar = os re Vo Gh) 
Vir) = Vo = Const. 


Pere or1sson's equation and the definition, of r.m.s. radius 


<r°> = fr-o(r)ar, (A-1) becomes 
Pee ore [> ones wee 
Deceviar = = o{(r)dr = “aan > (A-2) 


Equating the right hand sides of (A-1) and (A-2) yields, in 
leave onalized MKS units 
2 


ote 2. Un ake 


- 


ae 
where €) is the permitivity of free space and in (Fermi) 


Leg 





V (eV) x r- x Bne 
ao Oe OS 


Some 


2 


Ae oe 1072 Vv, (eV) ae (A-3) 


From (2-90) the expression for the charge form factor is 


given by 


aoe eo ae (A-4) 


pomenar the slope at the limit qé + Q gives 











dG, - 
- = om Zs) 2 oe x 107° Vo (eV) Bo Cheap 
dq q@+0 
Krohn and Ringo's [16] latest value is V, = -3720 eV, hence, 
dG, ; 
5 = +0.0193 F’. 
dq q©+0 
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APPENDIX B 


Moor nULio OF Tih TENSOR OPERATOR S45 


several experimental facts, like the nonadditivity of 


the magnetic moments (i.e., uw. + H # Ha) and the presence 


p 
of the quadrupole moment Qa in the deuteron, led to the 
mwuemesTtion that the nature of the nuclear force is 

eoeerally noncentral, that is, it contains a partly "tensor" 
hoOmce , 

Twitter cMay Fhe Censor operator be symmetric in 
mameeeres 1 and 2, and that it be invariant under the parity 
transformation (space reversal). Let r = 5a = Be be the 
frelauive coordinates of particle 1 and 2. The particles 


have spins given by o5 and o5 where o's are the Pauli spin 


matrices (an explicit representation) 


o = oOo = o_ = cae 


[o,, eye = e654 aod ra Ky VY 92 (B-2) 

oxo = 2id (B-3) 

of = of = of sank (B-4) 
x y Z, 


(6:h)(6-B) = (£-B) + ig: (A xB) (B-5) 


Peplyveation of the parity operator P to r and o yields, 


> > -> > 
Pr Te es PG a Ca 


deal 








Thus the vector Yr can appear only an even number of times 

in the tensor operator. Further, since the tensor operator 
: - -> -> 

iemeaymnetric in 1 and 2, Oy and O05 should appear in some 

bilinear combination. With (B-2,3,4,5) in mind it is found 

that the only linearly independent quantities which can be 


> > > 
made from r, 05 and 0, are 


> > 


> 
(oO, C5) and (o 


Be Coen) (B-6) 


The traditional definition then for the tensor operator 


Si is 


ompesece, © )(o, 7) = (o,-a,) (B-7) 


where r = r/|r|. 


The tensor operator (B-7) has some interesting proper- 


ties. namely <S4 5? Goyer all Cirecvrvons “sanesnes. Censader 
Ee A =e me Je => bees = JA ' 
<(o, r)(o, r)> = 2 fa, r) (0, r) sin@dé do (B-8) 


Take o4 along z-axis, o5 in the x-z plane making angle a 


ree rne 42 axis, and Y a vector making angles @, and @. 


Then 
eto (B-9) 
on = cetera x + cos a2) (B-~10) 
r= sin ¢@ cos $x + sin @ singy + cos @ 2 (B-11) 


Bienen (5-9,10,11) in (B-8) and integrating over @ and 9 


yields 


lia 


a 


a 





COS Qa = 


rae 
o~ 
Qy 


<(9,°r)(o, er = s 0,0. Geeks 


Hence it follows that <5) 5? = 0, 
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APPENDIX C 
CALCULATION OF DEUTERON MAGNETIC 
AND QUADRUPOLE MOMENTS 
The neutron and proton magnetic moments Wn? and Hy have 
been quite accurately determined, as has tetoee Cle ee teen 


magnetic moment Ha. The difference Stes + HW? _ = 0.0222rmm 


Md 
is small but not zero and much larger than the experimental 
perer. A possible explanation is that the deuteron spends 
a small part of the time in some state other than the 
symmetric ao Stave. The task is to find the correct mixture 
Meeteeres that Can account for the deuteron's moment. 

Recall the angular wave functions for the oa and =i 


States 


Se ieee Oral a 
a Pai cr ae 


oa i (So ee eo eo : 
Dy? 4 oni “V2 5 X1 io %2 X1 *Wio Yo Xz 62-57) 


Using (2-56) and (2-57) construct the deuteron wave function 


—p 
doaz(P). Tels 








Sen m Es 
ol eo Cae 
aoe. 1,2, _ wir) 1 _9 
Dy? $543 '") Saige OC 
The nuclear magnetic moment operator can be written 
Zi A 
tl = Q < e C= 
Mz > le : By S54) : En ok (C-3) 
k=1 k=Z+t+] 
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where dy is the orbital angular momentum of the ee pee LOM 


Sy Ghe spin angular momentum of the ee proton (or neutron) 
and where Ey and g, are Che oy romero rlommayios Of Lhe 
Proton and neutron, B, = Oe -3.83. (The carot over 
Ho denotes U, as an operator.) In the e-m system and for 


the deuteron (C-3) becomes 
-— 1 
aa ek St eS (Coatl) 


where the orbital angular momentum associated with the 
meocon 1s half of the relative orbital angular momentum. 
The magnetic moment is found by taking the expectation value 


Of (C-4) between the state dies 


= > *m «”“ m _ m A m 
<u? - fat boss Uy boss =<beszl¥gl oess? ae) 


The following formula are useful in evaluating (C-5) 


uYy 7 mY SD n Un on ones oelatg 
S 0X4 A i. Sonk] . : xy 

Sop x4 2 : Sonky = 7 . 3 (C-6) 
Mea 2 So ee 5 

Sap Xo wey no Son X1 = 7 : a 


re) 





Hence 


e me 1 il ome 
be Fort 2 (Er a 6.) 4 oi 5 2 (E, = 84/76 X11 


; iS oa ee ori 
be Foil Vi f- 2 )x AN SE fy Y2[x1tep-e,) : 


emo on y9 yl 
10 5 ae 


5 


(C-7) 


The orthonormality of the spherical harmonics Y integrated 


meee the solid angle, and the orthonormality of the spin- 


m 
wave functions ee applied to (C-7) yields 


36 E 


. 1) 1 
1: <b977 14,1937? 


On oo sein 


tl 


CO 


“ i 
i $577! 4 2! %17?= 3 7 i (e,t@,)] [wren 


0 


E T Ey eC t n 
a U (r)dr = a 
0 


(C-9) 


lt -F (e, + &,) | = oma cnet 


In evaluating the radial integrals the assumption was 


maaqe that the deuteron was in a definite state, either S or 


D but not both. But the deuteron magnetic moment is 
0.8573 nm, that is, neither of the two values above. 


®eonsider amixture of the two states. 


Thus 


iiemtonal wave function for a combination of states is 


1 
ene oD. tel te a 


5 
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CE=10)) 





Normalization gives 


<o|o> = ee (C-11) 
where ae =e u“(r)ar, nS Jf w°(r)ar are the individual 
0 0 


Probabilities that the deuteron is in an 5S or D state. Thus 


wes pee : 
"3543p ple, +E,iP.¢ [7 rE, +&,)]Pp (C-12) 


Pauations (11) and (12) are two simultaneous equations in 
the two unknowns i and Ph: Using the experimental value 


hOommlmirne (C—-l-) the solution for P. and P.. are 


oi D 


dq 


P= 4%, P_ = 96% 


S 


The quadrupole moment Qa can be found in the same manner 
as the magnetic moment, 


a Se Ae 


Qe = 37377 -r = (3 68° 


eee (C-13) 


In the c-m system, the distance of the neutron and proton 
Prom the center of mass is r/2, and only the proton con- 


teabutes to Qa: Hence, 


<Q,> = gO pOnae = apiece (c-14) 


In terms of spherical harmonics this becomes, 





Sih /167 ¥ 2.0 {> 
ei = 7 5 fc ie Yo aie 

iL Asean eee) 1 
cee a 161 [, “Yori 2 > Yor3+ oul a3 ee Jory 


¥ 1 2 oY 7 2 
Yon ots,2) + we Hoa t Yo Forte 
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3 
a 


integrated over the solid angle it is zero because of 


Mie tirst term represents the pure contribution and when 
spherical symmetry. The next terms in parentheses are 
Gernot an conjugates of each other and the total can be 


Written as 


; Sn vo ¥ ee 
uw 37 Dat aie 


which after angular integration gives 


1 | 2 
— / uvr 
¥50 ~5 


The third term comes from the 


3p 


1 state and gives 


co 


= Bf we? Gites 


0 


<Q? = Len’ = els ro (C-15) 
0 QO 


iemso further, explicit solutions for the radial wave func- 


Hence 


tions u(r) and w(r) are needed. Some examples of Q 4 for 


Ppectii1ce numerical solutions u and w are listed in Table III. 
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APPENDIX D 


THE BREIT OR BRICKWALL FRAME 


The Breit frame is a Lorentz frame in which either 


a = ft > >y . : , 
Ky As Ky or, Py + Py vanish, where in an elastic two-body 


t 


feilision, Kk,» Dy eee eciceni me cigh chews. —lOMem lila les 


D are outgoing particle 3-momentum. The somewhat odd 
memreion Stems from traditional techniques in elastic 
meavcering theory and dispersion relations, particularly 
the Mandelstam representation. Consider Figure D-l. The 
Pemvention is that all four-momentum are ingoing. This has 
the advantage that one may consider any two of the four to 
be the incoming particles and the other two as outgoing 
mies, the physical momentum of the outgoing particles is 


then the negative of the one used in Figure D-l. 


x Po 


Ky Py 


Figure D-1l 


The physical momentum of an outgoing particle is denoted 
by a prime. For example, if the particle corresponding to 
the arrow with Ky is outgoing, it has physical momentum 


1 


Ky =). 


ee 





Mopeyea Lorencz transformation to the scattering event 
such that K, 1° kK. = 0. Thus, four-momentums Ky and kK have 


m@ewrorm (h = c= 1) 


k, = (wk), ky = (w,-k) (D-1) 


From conservation of energy it follows that the energies of 


Meemomparticle before and after the collision, Ey and 2a 


must be equal. Hence Ip, | = Ip. | and 
{ 
BEES eP,)> Po = (€,p,) Wace 
with 
IB,| = [Bol =p = yee - 
Da P5| =p 


Since scattering is elastic k, - k,, =p, - p, and it follows 


1 


k, -~ ki} = (0,2k) = ph - py = (0,B5-D,) Giese) 


2k the three momentum transfer. Hote wonmmel) (2) vand (3) 
mela the following picture. Both particles seem to be 
reflected on a hard wall, the k particle perpendicularly and 


the particle at some angle 





Figure D-2 


re 





In the c-m system the energies € and w of the p and k 
particles are not independent, but they are in this system. 
Therefore they are convenient to use as variables. In 


this system the four momentum transfer is very simple: 


2 


q = (k, - k3)° KR) 


>) = (0,2k)° = - |2k|° =-q 


mew, 1b 1S the three momentum transfer. 
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